Approximate Waveforms for Extreme-Mass-Ratio Inspirals 
in Modified Gravity Spacetimes 
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Extreme mass-ratio inspirals, in which a stellar-mass compact object spirals into a supermassive 
black hole, are prime candidates for detection with space-borne milliHertz gravitational wave detec- 
tors, similar to the Laser Interferometer Space Antenna. The gravitational waves generated during 
such inspirals encode information about the background in which the small object is moving, provid- 
ing a tracer of the spacetime geometry and a probe of strong-field physics. In this paper, we construct 
approximate, "analytic-kludge" waveforms for such inspirals with parameterized post-Einsteinian 
corrections that allow for generic, model-independent deformations of the supermassive black hole 
background away from the Kerr metric. These approximate waveforms include all of the qualita- 
tive features of true waveforms for generic inspirals, including orbital eccentricity and relativistic 
precession. The deformations of the Kerr metric are modeled using a recently proposed, modified 
gravity bumpy metric, which parametrically deforms the Kerr spacetime while ensuring that three 
approximate constants of the motion remain for geodesic orbits: a conserved energy, azimuthal an- 
gular momentum and Carter constant. The deformations represent modified gravity effects and have 
been analytically mapped to several modified gravity black hole solutions in four dimensions. In 
the analytic kludge waveforms, the conservative motion is modeled by a post-Newtonian expansion 
of the geodesic equations in the deformed spacetimes, which in turn induce modifications to the 
radiation-reaction force. These analytic-kludge waveforms serve as a first step toward complete and 
model-independent tests of General Relativity with extreme mass-ratio inspirals. 

PACS numbers: 04.30.-w,04.50.Kd,04.25.-g,04.25.Nx 
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I. INTRODUCTION 



"I seem to have been only like a boy playing on the 
seashore, and diverting myself in now and then find- 
ing a smoother pebble or a prettier shell than ordinary, 
whilst the great ocean of truth lay all undiscovered be- 
fore me." [l[ Isaac Newton's quote reminds us of the great 
unknowns that remain in gravitational astrophysics Q. 
Our understanding of this field should soon be revolution- 
ized by the detection of gravitational waves (GWs) with 
ground [!, H| and space- based interferometers (H-IU . In- 
deed, the ground-based detectors LIGO [1] and Virgo [To[ 
are currently undergoing upgrades toward a sensitivity 
at which a first direct detection of gravitational waves is 
likely. Space-borne detectors are being planned and will 
hopefully be operational in the next decade. 

Detectors in space will be sensitive to low-frequency 
[(10 _5 -10 _1 ) Hz] GWs. One of the most promising 
sources in this frequency range are extreme-mass ratio in- 
spirals (EMRIS) [ll]]. These systems consist of a stellar- 
mass compact object [(10°-10 2 )M Q ], such as a neutron 
star or black hole (BH) , that spirals into a supermassive 
BH [(10 5 -10 7 )Af Q ], typically on an inclined and eccentric 
orbit. Due to their extreme mass-ratio, these inspirals 
proceed slowly, generating hundreds of thousands of cy- 
cles of gravitational radiation while the smaller object is 
in the strong-field region close to the central supermas- 
sive black hole. These GWs encode detailed information 
about the structure of the spacetime exterior to massive 
compact objects and the non-linear, "strong-field" nature 



of the gravitational theory that describes the dynamics 
of the inspiral. 

The detection of such GWs requires the construction 
of accurate waveform templates that allow the extrac- 
tion of signals from noisy data via matched filtering. 
EMRI waveforms, however, are particularly difficult to 
construct, as one requires these to be accurate over hun- 
dreds of thousands of cycles. Progress has been made 
by constructing approximate waveforms that encode all 
the key features of EMRI waveforms, while being com- 
putationally inexpensive. There are two families of such 
"kludge" waveforms that have been used extensively — 
the "analytic kludge" developed by Barack and Cut- 
ler (l2j . which we will focus on in this paper, and the 
numerical kludge [H|,[l4[- In the numerical kludge, the 
orbital trajectories are based on exact Kerr geodesies, the 
parameters of which are slowly evolving under the influ- 
ence of radiation-reaction. Once the orbital trajectories 
have been obtained, the waveforms are constructed us- 
ing a quadrupole approximation. The AK model is built 
around the gravitational waveforms generated by parti- 
cles describing Keplerian ellipses [HI, ED], whose semi- 
major axis, eccentricity and inclination angle evolve ac- 
cording to certain post-Newtonian (PN) equations (l~5l — 
[l7| . Relativistic precession of the orbital plane and the 
perihelion are also included using post-Newtonian ap- 
proximations. 

In the past, attempts have been made to modify EMRI 
waveforms to study their ability to test General Relativ- 
ity (GR). Such attempts can be divided into two classes: 
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"extrinsic" or "intrinsic" . The former picks a concrete al- 
ternative theory, such as dynamical Chern-Simons (CS) 
gravity jla. ll9l , constructs waveforms for that particular 
theory [20l - |23j and compares these to the predictions of 
GR. The latter are null-tests of GR, where one assumes 
the validity of GR a priori and concentrates on internal 
consistency tests, such as measuring the multipolar struc- 
ture of the metric [23 - |27j , or multi- modal spectroscopy 
of inspiral and ringdown waveforms (28l - [30| . Both classes 
of tests are intrinsically valuable, but they are not ideal 
to search for departures from GR in a systematic and 
model-independent way. 

Recently, there has been a focused effort to develop 
such a systematic and model-independent approach. 
Yunes and Pretorius [3ll - l33j proposed the parameterized 
post-Einsteinian (ppE) framework, a model-independent 
template family for the quasi-circular inspiral of com- 
parable mass, non-spinning, compact binaries. In this 
approach, deformations of the conservative Hamiltonian 
and of the radiation-reaction force are mapped onto the 
waveform observable, i.e. the frequency-domain GW re- 
sponse function. The meta-template family allows for 
generic deformations away from GR, that have been 
shown to reproduce waveform predictions from all known 
alternative theories proposed to date. However, this ap- 
proach is ill-suited to EMRI waveforms, as the orbits are 
likely to be inclined and eccentric and the mass ratios are 
extreme. 

To address this shortcoming, Vigcland, Yunes and 
Stein [34| took the first steps toward extending the ppE 
framework to EMRIs. They concentrated on conserva- 
tive corrections to the orbit, realizing that these could 
be parameterized by deformations of the metric ten- 
sor. Previous attempts to construct such "bumpy space- 
times" [25U27I had focused on intrinsic or null tests, 
and had therefore used metrics that satisfied the Ein- 
stein equations, but with potentially unphysical matter 
distributions, such as naked singularities, or unphysical 
spacetime regions, e.g., closed timelike curves. The ap- 
proach of [34[ was to only allow for the subclass of met- 
ric deviations that would ensure the existence of an ap- 
proximately conserved energy, (z-component of) angular 
momentum, a second-order Killing tensor and a Carter 
constant, without requiring that the Einstein equations 
be satisfied. Such an approach led to a paramctrically 
deformed metric that was shown to map to metrics in 
known alternative theories in four dimensions [20 . |35| . 

In this paper, we construct corrections to the AK wave- 
forms by considering geodesies in the modified gravity, 
bumpy metric of 3J]. We begin by providing explicit 
expressions for all components of the metric deformation 
as an expansion in r 3> M, where r is the field-point 
distance to the supermassive BH with mass M; we re- 
quire that the metric remain asymptotically flat, with 
the same scaling at spatial infinity as that predicted by 
the GR peeling theorems, and disallowing pure angular 
deformations. We show that even with these restrictions, 
all known metrics in alternative theories arc still recov- 



ered by an appropriate choice of the parameters charac- 
terising the metric functions [2(| H^] . Arbitrary choices 
of these parameters lead to parametric, metric deforma- 
tions at C(l/r 2 ), 0(l/r 3 ), 0(l/r 4 ) and C(l/r 5 ) relative 
to the Kerr metric. 

We then calculate the geodesies equations associated 
with this background. We parameterise the orbits via 
the location of their turning points and will relate vari- 
ous quantities to the Kerr geodesic orbit with the same 
turning points. The existence of three approximately 
conserved quantities allows us to explicitly separate the 
geodesic equations into first order form. From these, we 
calculate the three orbital frequencies associated with 
the orbital motion. Although we do not introduce ex- 
plicit corrections to the radiation-reaction force, modifi- 
cations to the fluxes of energy, angular momentum and 
Carter constant will be introduced due to modifications 
to the orbit itself. We calculate these implicit modifi- 
cations with a quadrupolar approximation for the fluxes 
of energy and angular momentum. The Carter constant 
flux is calculated by assuming that the inclination angle 
remains approximately constant and by requiring that 
exactly circular orbits remain circular under radiation 
reaction. Finally, we collect all pieces of the calculation, 
providing an explicit prescription to build AK waveforms 
in these families of parametrically deformed spacctimes. 

The study performed here is similar, yet more generic 
than others already carried out in the literature. For 
example, Barack and Cutler [28| considered modifica- 
tions to AK waveforms induced by a perturbation to the 
quadrupole moment of a Kerr BH. Using results in [36j , 
they introduced modifications to the precession frequen- 
cies and rate of change of orbital frequency that would 
be induced by a quadrupole moment deformation; they 
then searched for the accuracy with which LISA could 
measure the size of such a deformation. This study, how- 
ever, neglected modifications to the eccentricity evolu- 
tion. Glampcdakis and Babak [37| also considered a class 
of bumpy spacetimes which differed from Kerr in the 
quadrupole moment only. Unfortunately these metrics 
are of Petrov Type I, and thus do not allow for the exis- 
tence of a Carter constant or for the separability of the 
geodesic equations |38j . Moreover, both studies neglected 
the possibility that the metric could be modified at mul- 
tipolc orders higher than quadrupole. This is a criti- 
cal disadvantage, as strong-field modifications of GR arc 
likely to introduce corrections at higher than quadrupole 
order, e.g., dynamical CS gravity modifies the metric at 
hexadecapole order, leaving the quadrupole and octopole 
unchanged. 

Although the waveforms presented here will be useful 
for studies that determine the accuracy to which instru- 
ments like LISA could constrain model-independent de- 
viations from GR, they are not unique or complete. To 
demand uniqueness is futile in a ppE scheme, as it is also 
in the parameterized post-Newtonian or parameterized 
post-Kcplcrian schemes that are used to search for devi- 
ations from GR in the Solar System and in binary pulsar 
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observations respectively. An infinite number of theories 
predict an infinite number of possible deviations in the 
observables of interest, while the schemes above can only 
parameterize a subset of them, i.e., the subset that can 
reproduce all the predictions of alternative theories that 
are known to date [2(], IH] . 

These waveforms are also not complete because here we 
will modify only one of the three main ingredients that go 
into waveform construction, i.e., the metric tensor, and 
will neglect modifications to the first-order equations of 
motion (that prescribe how GWs are sourced by mat- 
ter distributions) and the second-order equations of mo- 
tion (that describe the self- force and radiation- reaction). 
However, a large sub-class of quadratic gravity theories 
exist in which the corrections to the metric are dominant 
over modifications to the wave generation and radiation- 
reaction J3j|. Nonetheless, a more complete analysis 
should also investigate the excitation of scalar and vecto- 
rial modes in the metric perturbation, which could arise 
from modifications to the wave generation. We leave an 
investigation of these other effects to future work, but 
the results found here will still hold when these other 
modifications are introduced. 

The remainder of this paper is organized as follows: 
Section [TT] introduces the new bumpy framework to 
model modified gravity theories; Section IIIII calculates 
the geodesic equations that preserve the Kerr turning 
points; Section IIVI computes the orbital frequencies as- 
sociated with the modified orbital motion; Section [V] 
calculates the implicit deformations introduced to the 
radiation-reaction force; Section lvTl builds AK waveforms 
from the modified orbital frequencies and fluxes; Sec- 
tion IVlTl concludes and points to future research. 

Throughout this paper we use the conventions of Mis- 
ner, Thorne and Wheeler [Z(|. We use Greek letters to 
denote spacetime indices, while Latin ones in the middle 
of the alphabet stand for spatial indices only. We 

also use geometric units with G = c = 1. Background 
quantities are denoted with an overhead bar, while quan- 
tities associated with geodesies in the Kerr metric are 
denoted with a subscript K. 



II. BUMPY SPACETIMES FOR MODIFIED 
GRAVITY 



In this section we discuss the parametrically deformed 
BH metric that we will later use to study geodesies. We 
begin by recapitulating the most important results of (34j 
for this paper. We then simplify this metric prescription 
by considering expansions in M/r <C 1, thus allowing us 
to compute explicit expressions for all components of the 
metric deformation. Finally, we describe some properties 
of the new metric. 



A. Spacetime Construction 

We decompose the metric tensor that is to describe the 
background spacetime of a supermassive BH as 



(1) 



where e < 1 is a "deformation" book-keeping parameter 
that reminds us that l^i/l/lg^i/l <C 1. The background 
metric is assumed to be the Kerr metric g^ v = , which 
in Boyer-Lindquist coordinates has components: 
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for a BH with mass M and spin angular momentum di- 
rected along the symmetry axis of magnitude S = M 2 a, 
where a is the dimensionless Kerr spin parameter. Equa- 
tions ©-© depend on the functions 



p 2 = r 2 + a 2 M 2 cos 2 9 , 
A = r 2 f + Af 2 a 2 , / = 1 
£ = (r 2 + M 2 a 2 ) 2 - M 2 a 2 A sin 



2M 



(4) 
(5) 
(G) 



We restrict attention to a certain class of metrics. 
We begin by requiring that the full metric be station- 
ary and axisymmetric, although it need not solve the 
Einstein equations, i.e. it is a solution to a more gen- 
eral set of modified gravity field equations that have a 
smooth GR limit: g^ v —¥ g^u a s hau 0. In addi- 
tion to the existence of a temporal and an azimuthal 
Killing vector, we also assume that a certain integrabil- 
ity condition holds (see Eq. (49) in [iH) such that the 
metric can be written in Lewis-Papapctrou form. The 
deformation of this metric g^ = g^ + e h^, can be 
transformed to Boyer-Lindquist-like coordinates in which 
all components of the metric perturbation vanish except 
{htt,ht$,h rr ,h r Q,hQ8,h^^). These are the only compo- 
nents that are allowed to be non-zero. 

With this at hand, we then force the full metric g^ v 
to possess three constants of the motion: a conserved 
energy, azimuthal component of angular momentum and 
Carter constant. The first two are generated directly 
from the Killing vectors and are exact, while the last one 
is built from an approximate, second-order Killing tensor 
valid at least to 0(e 2 ). This Killing tensor is parameter- 
ized as 



(7) 



The parameterization of the Killing tensor as in Eq. ([7]) 
identifies the coordinate r with a Boyer-Lindquist-like ra- 
dial coordinate. The vectors and are required to 
be null and is required to satisfy the Killing tensor 
equation, which determines the components of the null 
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vectors up to some arbitrary functions of the radial co- 
ordinate. This in turn determines the final form for the 
metric perturbation. 



The non- vanishing components of the metric deforma- 
tion can be written as [3JI Egl 
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where p 2 = r 2 — M 2 a 2 cos 2 and PP K are polynomials in 
r and cos 0, given explicitly in the Appendix of [34[ (we 
have here adopted the deformed Kerr parameterization 
of [34j|). The quantities 7, = 7<(r) are arbitrary functions 
of radius, while 03 = 63(0) is an arbitrary function of 
polar angle. The quantity I is defined as 



1= dr 



2M 2 a 2 sin cos 

2 ' 

p 2 



L r9 



2r 71 



dr 



(13) 



The metric perturbation component h r8 is free. 

The parametrically deformed metric represents a fam- 
ily, some members of which are well-known BH solu- 
tions in modified gravity theories. For example, a certain 
choice of deformation parameters 7^4 leads to the slowly- 
rotating BH solution in dynamical CS gravity [2(| , while 
another choice leads to modified Schwarzschild BHs in 



quadratic gravity theories [35||, as was shown in [34| . In 
both cases, these solutions derive from field equations 
that arise from a diffeomorphism invariant theory, with a 
well-defined Lagrangian density. We refer the interesting 
reader to 1341 for more details on the metric construction. 



B. Simplification of the Parameterization 

Let us simplify the metric perturbation of the previous 
section with the following criteria: 

1. Asymptotic flatness: Require that h^ v — > at 
spatial infinity. 

2. Peeling: Require that \hy, v \ ~ r -2 or faster for 
r/M > 1. 
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3. Occam's Razor: Set the largest number of metric 
components to zero that are not needed to repro- 
duced known modified gravity predictions for the 
metric tensor in four dimensions [20l . |35| . 

Requirement (1) ensures that there is no constant piece 
to the metric deformation, such that the total metric 
9\-iv ~~ ^ 9nv a t iP ■ Since the Kerr metric is asymptoti- 
cally fiat, so would the total metric. Requirement (2) 
ensures that the BH mass is not renormalized by 1/r 
corrections to the (t, t) or (r, r) metric components. The 
norm in this requirement is to be taken with the flat 
metric in spherical coordinates, such that the {9, <f>) sub- 
sector is simply the metric on the 2-sphere. These two 
requirements ensure that the metric perturbation satis- 
fies all Solar System constraints, as it introduces mod- 
ifications at higher than leading, Newtonian-order in a 
weak- field expansion. Furthermore, any 1/r correction 
to (htt, httj), h rr ) would renormalize the BH mass or spin 
angular momentum, which would not be measurable. Re- 
quirement (3) automatically implies that (i) we can set 
h r = 0; and (ii) we can disallow pure angular deforma- 
tions by setting 63 = 0. 

Since the purpose of this paper is to construct ppE 
AK waveforms, which make extensive use of weak-field 
expansions in M/r < 1, we choose to parameterize the 
remaining free functions as Taylor series: 




(14) 

where A = 1 or 4 and 7^ n are dimensionlcss constants. 
We have pulled out a factor of 1/r in the expansion of 
73 because this quantity has dimensions of [M] _1 , as one 
can see from in Eqs. (|5j)- (fT2"j) . Notice that with these 
definitions the 7 TO , n constants are dimensionless. 

With this at hand, let us simplify the metric compo- 
nents, starting with h rr . The integral I can now be 
solved exactly: X= p 71. Via requirement (3) 

P 2 

h rr = --y 1 — . (15) 



Requirements (1) and (2) force us to choose 71.0 = = 
71,1, since p 2 /A — > 1 for r 3> M. 

The next simplest component to analyze is hgg, whose 
behavior is governed by Eqs. (|5|)- ([T2")l , which using the 
previous results simplifies to dh gg /dr = 2r / (p 2 )h gg . 
The solution to this equation is hgg = <d±(9)p 2 . Since 
999 = p 2 , this correction would be leading order in the 
angular sector. By requirement (3), we disallow it and set 
04 = 0. This simplification, and the restrictions made 
when deriving the metric perturbation, fix the coordi- 
nate system in such a way that hgg = and therefore 
9ee = p 2 +0(e 2 ). The radial coordinate thus preserves, at 
leading order, its physical interpretation in the Kerr met- 
ric — it may be interpreted as a circumferential radius in 
the equatorial plane and becomes the oblate-spheroidal 
radial coordinate in flat-space at infinity. The fact that 
the radial coordinate has been fixed in this way will be 
important for interpretation of the waveforms we derive 
in this paper. This will be discussed further in ScctionfVT] 

To determine the remaining metric components, we 
must first solve the differential equation for h t $, as the 
htt and h^ components depend explicitly on the former. 
This is an elliptic equation that could be solved numeri- 
cally. Since we seek analytic solutions only, however, we 
will solve it in the M/r <C 1 limit. To do so, we write 




where ht^n are functions of 6 that we will determine by 
solving Eqs. (|5])- ([T2"|) . We could have chosen to expand 
httj, in a 1/r basis, but we empirically found that a 1/p 
basis yields simpler results. Solving Eqs. ([5T)- (fT2")) order 
by order in 1/p, we find that the first few non-zero terms 
are 



h t <t>,2 = sin 2 6 [-73,3 - a (71,2 + 74,2) - a 2 73,i] > ( 17 ) 
= sin 2 9 [(273,3 - 73,4) + a (674,2 - 74,3 + 271,2 - 7i.s) + 2a 2 73,i (4 cos 2 9 - 3)] (18) 

h t <p,4 = sin 2 9 [-a 4 73,i - (74,2 + 71,2) a 3 + 2 (-473,1 cos 2 9 - 73,3 + 473,1) a 2 

+ (-71,4 - 874,2 + 674,3 - 74,4 + 271,3) a + 273,4 - 73,5] , (19) 

ht*,s = sin 2 9 {273,5 + a (-74,5 - 874,3 + 2 7m + 674,4 - 7i,s) + a 2 [-273,4 + 73,3 (-2 + 5 cos 2 9) - (1/2)73,4 cos 2 9} 

+ a 3 [(1/2) cos 2 9 (271,2 - 74,3 - 7i.3 - 274.2) - 74,3 - 7i,3 + 474.2] + a 4 (-473.1 cos 4 9 + 973,1 cos 2 9 - 473,1) } , 

(20) 

ht^fi = sin 4 9 [a 2 (73,5 - 673,4 + 873,3) + a 3 (71,4 - 274,3 + 74,4 - 271,3) + 273,3a 4 

+ a 5 (74,2 + 71,2) + a 6 73,i] + sin 2 9 [273, 6 - 737 + a (674,5 + 271,5 - 874,4 - 7i,6 - 74,6) 
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+ a 2 (473,4 - 373,5) + a 3 (-271,4 + 271,3 + 674,3 - 274,4) - 373,3a 4 - a 5 (71,2 + 74,2) - 73,i a6 ] > ( 21 ) 



where we have here simplified the solution by setting 
74,0 = 0, 74,1 = and 73,2 = 0. We will find that these 
conditions are necessary to ensure the metric is asymp- 
totically flat. We have also set 73,0 = 0, which is required 
for the differential equation to be satisfied. 

Let us now return to the h tt and components, 
which have been completely specified by the above so- 
lutions. At spatial infinity we find that 



2 7 , 



4,0 



274,1 - 8M 7 4 i0 



.^73 l2 sin^ + 0(l) 



(22) 
(23) 



when M/r <C 1. By requirements (1) and (2), this then 
implies that (74,0,74,1,73,2) must all be set to zero. 

In summary, the requirements of asymptotic flatness 
and the non-renormalization of the mass, have forced us 



to the following conditions 
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, 71,0 
, 74,0 



0, 
0, 



71,1 
74,1 



0, 
0, 



73,0 
73,2 



(24) 
(25) 



With these choices, the metric perturbations (hgg,h r g) 
vanish, h rr is given by Eq. (1151) and h t< j, is given in the 
far field by Eq. (|16[) , with the angular functions given in 
Eq. (|21[) . The remaining components (htt-, h^) are given 
explicitly by Eqs. (|g|)- ([T2)) . 

Let us now take the far-field expansion of all the metric 
perturbations: 



(26) 



The first few non-zero terms are 



htt, 2 = 71,2 + 274,2 - 2a73,i sin 2 9 , 

htt,3 = 7i,3 - 874,2 - 271,2 + 274,3 + 8073,1 sin 2 6 , 

htt a = -874,3 - 271,3 + 274,4 + 874,2 + 71,4 - 8073,1 sin 2 9 + a 2 (71,2 + 274,2) sin 2 9 + 2a 3 7 3 ,i cos 2 6 sin 2 9 , 
htt,s = 16a 3 7 3 ,i sin 4 9 + sin 2 9 [4a 7 3, 3 + a 2 (71,3 - 2 7l)2 - 12 7 4, 2 + 2 74 , 3 ) - 12a 3 73 ,i] 

+ a 2 (874,2 + 271,2) + 71,5 + 274,5 - 271,4 + 874,3 - 874,4 , (27) 



h rr ,2 = -71,2 , h rr , 3 = -71,3 - 271,2 , h rrA = -71,4 - 271,3 - 471,2 + (1/2)71, 2 a 2 (1 - cos26») , 

h rr ,5 = a 2 sin 2 6 (71,3 + 271,2) - 71,5 - 271,4 - 471,3 - 871,2 + 2a 2 7i, 2 , (28) 

h t <t>,2 = —M sin 2 9 [73,3 + a (71,2 + 74,2) + a 2 73,i] , 

ht<f,,3 = -8Ma 2 73,i sin 4 9 + M sin 2 6 [(273,3 - 73,4) + a (674,2 - 74,3 + 271,2 - 71.3) + 273,1a 2 ] , 
h<f,A = Msin 4 9 [a 2 (873,1 - 73,3) + a 3 (-71,3 - 74,2) - a 4 73,i] + sin 2 9 [(273,4 - 73,5) 

+ a (-74,4 - 874,2 + 674,3 - 7i,4 + 271,3) - a 2 73, 3 ] 
htfrS = -16Afa 4 7 3,i sin 6 6 + A I sin 4 9 [a 2 (-273,3 - 73,4) + a 3 (-74,3 + 1074,2 + 2 7 i, 2 - 71,3) + 14a 4 7 3 ,i] 

+ sin 2 9 [(273,5 - 73, 6 ) + a (-71,5 - 874,3 - 74,5 + 271,4 + 674,4) - 73,4a 2 

+ a 3 (-271,2 - 674,2) - 2a 4 73 ,i] , (29) 

h^-2 = , h^-i = , = 2M 2 a7 3 ,i sin 4 6 , = , 

h^, 2 - M 2 sin 4 6 [2a 7 3, 3 + a 2 7i,2 + a 3 73,i (4-2 cos 2 0)] , 

hwfl = 8M 2 a 3 7 3,i sin 6 + M 2 sin 4 [a (-4 73 ,3 + 2 73 , 4 ) + a 2 (-2 7 i, 2 - 4 74 , 2 + 71,3) - 4a 3 73 ,i] , (30) 

I 



Notice that although (h t t,h rr ) are indeed dimension- squared, as expected since these are also the dimensions 
less, httf, has units of length and h^ has units of length 
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of the corresponding components of the Kerr metric. 

The perturbation is parameterized by a number of 
constants, depending on how many terms in M/r are 
kept relative to the leading-order Kerr metric: up to 
0(M 2 /r 2 ), the metric deformation is given by the 4 
constants B 2 = (71,2, 73, 1, 73,3, 74,2); up to 0(M 3 /r 3 ) 
it is given by the 7 constants B 2 U S3, where B3 = 
(71,3,73,4,74,3); up to 0{M 4 /r 4 ) it is given by the 10 
constants B 2 VJBz\JBi, where 84 = (71,4,73.5,74.4); up to 
0(M 5 /r 5 ) it is given by the 13 constants B 2 UB 3 UB 4 UB 5 , 
where B$ = (71,5,74,5,73,6)- Later on in the paper, we 
will take certain Bn limits, where we mean we will only 
let the Bn coefficients be non-zero and set all others to 
zero. 

Lastly, note that known BH solutions in alternative 
theory of gravity can be reproduced with this parame- 
terization, as shown in [34|. In particular, in this paper 
we will frequently compare our results to that of dynam- 
ical CS gravity, where the metric of a slowly rotating BH 
is identical to Kerr, except in its (i, <f>) component, which 
is given by [2(| 



9t<t> 



2M 2 ar 



A 



aM 5 



sm 



1 



12M 27M 2 



7r 



10r 2 



sin 2 e , 



(31) 



where £ is a dimensionlcss coupling constant of the the- 
ory. This metric can be completely reproduced by the 
above parameterization by choosing all 7^ = 0, except 
for 73 (r), whose first non- vanishing terms are 



73,5 



-a( , 73,6 



65 
28 



«C ■ 73,7 



709 
112 



oC- (32) 



C. Physical Properties of the Deformed Metric 

At this junction, one might wonder how the metric of 
Eqs. (|8"])- (fl"2|) may change physical properties of the space- 
time. For a stationary and asymptotically flat spacctime, 
the event horizon coincides with the Killing horizon, and 
it is given by the hypersurface where: 



at its Kerr value. Moreover, one can also show that 
the horizon coincides with the location of a coordinate 
singularity, by evaluating the components of the inverse 
metric. 

Other quantities, however, will be different in the de- 
formed spacetime relative to their Kerr values. For ex- 
ample, it is likely that the geometry of the ergosphere, 
defined by the hypersurface where gtt = 0, will be modi- 
fied. These and other properties of the spacetime metric 
will not be needed in this paper, and so we leave them 
to be explored elsewhere. 



III. GEODESIC EQUATIONS 

In this section, we study the geodesic equations asso- 
ciated with the new metric we computed in the previous 
section. We begin by summarizing the most important 
results of [34|, regarding the geodesic equations. We then 
fix the location of the turning points of the orbit and com- 
pute the modifications to the conserved quantities rela- 
tive to the Kerr geodesic with the same turning points. 



A. First-Order Equations 

Let us define the dimensionless constants of the motion 
associated with this new metric as follows: 



E 



-t a u a , L = M(j) a u a C = M 2 i af} u a u p , (36) 



where t a = [1,0,0,0] is a timclike Killing vector, <f> a = 
[0,0,0,1] is an azimuthal Killing vector and u a — 
[i,r,6,4>] is the 4-velocity, with overhead dots stand- 
ing for differentiation with respect to proper time. The 
Carter constant is also often defined as Q = C — (L — 
aE) 2 . Here the approximate, second-order Killing ten- 
sor £ a /3 = £ a /3 + S^ a p . The background Killing tensor is 
simply the Kerr one 

£ al 3 = Ak(Jf}) +r 2 g afS , (37) 

with k a and l a the principal null directions of Kerr: 



„2 



9u 



0. 



(r 2 +M 2 a 2 ) 2 

hrhrh . ,» „ lltt 



2(r 2 + M 2 a 2 



(33) 



Let us now use the fact that at r = r + = M + M (1 
a 2 ) 1/2 , A = and then 



ht* ■ (34) 



M 2 a 2 a 
With this at hand, linearizing Eq. in e leads to 







-A 



2ht 4 



Ma 



M 2 a 2 



(35) 



This quantity clearly vanishes at r = r + , which then im- 
plies that the event horizon of the full metric remains 



k a = 

r = 



r 2 + M 2 a 2 Ma 

A ' 
r 2 + M 2 a 2 Ma 
^A— '-^^ 



The Killing tensor deformation is 

6£ a p = A 8k( a lp) + 8l( a kp) + 2h S{a kp)l i 
where the deformed null vectors are 

„2 , J#2„2 



8k a 



Sl a = [74,0,0,73] 



(38) 
(39) 



3r 2 h 



r~ + M 2 a l aM 
^ 7i +74,71,0,^-71 +^3 



af) ■ 

(40) 

(41) 
(42) 
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and we have used the requirements of Sec. lIIBl to simplify 
these expressions. 

The dimcnsionless constants of the motion, when eval- 
uated on the perturbed metric, become 

E = E + e (h^tPu" + g^Su") , (43) 
L = L + Me {h^u v + g^Sv?) , (44) 

C = C + M 2 e (s^ v u»u u + 2^uHu v ^ , (45) 
and the normalization condition for the four-velocity is 



= KuU^u" + 2g ilv vf 1 Su v 



(46) 



since by definition — 1 = g flv u fl u 1 ' '. In all of these equa- 
tions, « M = [t, f, 0, 4>] is the unperturbed, Kerr four- 
velocity, while 8u^ is a perturbation of 0(e). 

The second-order geodesic equations can be rewritten 
as a first-order set through Eqs. (|4*3")l -(|4*5 ]) . To achieve 
this, one must make a gauge choice for the constants of 
the motion associated with the full spacetime (E,L,C). 
In [34[, the authors chose to keep E = E, L = L and 
C = C, which then implies that 8u^ must be such that 
all terms in parentheses in Eqs. (j43 ]) - (|45|l vanish. This 
choice, however, forces the turning points of the orbit to 
be different from the GR orbit with constants (E, L, Q). 

Using this condition and Eq. (|4l)]) . the geodesic equa- 
tions can be rewritten in first-order form: 



p 2 i = T K (r, 0) + eST{r,0) 
R K (r) + eSR(r,0), 



4 -2 

p r 



p 2 <j> = <J> K (r, 0) + e <J$(r,i 



(47) 
(48) 
(49) 
(50) 



p-r =e K (0) + e8Q(r,0), 
where the Kerr potentials (T K , R K , K , 'I'k) are given by 

P 



-aM (aME sin 2 - ML) + (r 2 + M 2 a 2 ) 



A ' 
(51) 



Rk = P 2 



M 2 Q + M 2 (aE-L) 2 +r 2 



(52) 



e K = M 2 Q - M 2 L 2 cot 2 - M 2 a 2 cos 2 (l - E 2 ) , 
ML \ aMP 



$ K = - aME - 



sin 2 



A 



(53) 
(54) 



where P = E(r 2 + M 2 a 2 ) - M 2 aL. The perturbation 
potentials (ST, 5R, SQ, (5$) are given by 



ST(r,0) 



(r 2 + M 2 a 2 ) 2 



M 2 a 2 sin 2 I 



h tn u a 



2aM 2 r 
A 



h ~ a 



SR(r, 0) = A [A(r, 0) r 2 + B(r, 0)} , 
SQ(r,0) 
6<5>{r, 0) 



(55) 

(56) 
(57) 



A(r, 0) M z a z cos 2 - B(r, 0) , 
2aM*r _ a p 2 - 2Mr _ a 

htaU 1 . 2 „ h <t><* U 1 ( 58 ) 



A sin^ 



and the functions A(r,0) and B(r,0) are given by 



A(r, 0) = 2 



h a tt + h a< f 



B{r, 0) = 2 Zt*$) Su* + (ttJ+M) 



(59) 
(60) 



in which the four- velocity associated with the background 
trajectory m m is 



u 
u 6 



t = p- 2 T K (r,0) u r = f = p- 2 ^R K {r), (61) 



P~ V@k(0) . u^ = ^ = p-H K (r,0), (62) 



and the four-velocity associated with the perturbation is 



Su* = St = p- z 5T , Su* = 5<f> = p~ z S<$> . (63) 



Let us expand the perturbation to the geodesic equa- 
tions in the far-field limit, using the metric in Eqs. (|27[1 - 
Doing so, we find 



00 / M\ n 00 / M\ n 

n=0 ^ ' n=2 ^ ' 



SR = M 4 J2 5R " 



n=-2 



M 



(64) 

8Q = M 2 jr SQ n 

(65) 



n=0 



where we note that (ST n , SR n , 5Q n , S<fr n ) are all dimcn- 
sionless, unlike (T K , ST, K , 8<d) which have units of M 2 , 
(R K ,SR) which have units of M 4 and ($ K ,<5$) which 
have units of M. With this, the first non-vanishing per- 
turbations arc 



ST = (274,2 + 71,2 - 2a7 3j 



1 sm 



E, 



STi = (71,3 + 274,3 + 271,2) E . 



8T 2 = E [(71,4 + 271,3 + 274,4 + 471,2) + a 2 (71,2 + 274,2)] + L [-73,3 ~ a (74,2 + 71,2) - a 2 73, 
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ST 3 = E [(71,5 + 274,5 + 471,3 + 271,4 + 871,2) + (71,3 + 274,3) a 2 ] - L [73,4 + (71,3 + 271,2 + 74,3) a] , 
ST 4 = E [(-471,2 + 274,4 + 71,4) a 2 + (274,6 + 871,3 + 271,5 + 1671,2 + 71,6 + 471,4)] 

+ L [-a 2 73, 3 - 73,5 - (74,4 + 71,4 + 471,2 + 271,3) a] , 
ST 5 = E [(-471,3 + 274,5 - 1671,2 + 71,5) a 2 + (47i, 5 + 3271,2 + 871,4 + 274,7 + 71,7 + 27i, 6 + 1671,3)] 

+ L [-73,4a 2 - 73,6 + (-471,3 - 74,5 - 71,5 - 271,4 - 871,2) a + 2a 3 7i, 2 ] , (66) 

5R-2 = (274,2 + 271,2) E 2 - 2^X73,! - 71,2 , 
6R-1 = E 2 (-474,2 + 274,3 + 271,3) + 4SL73,i + 271,2 - 71,3 , 
SRo = E 2 [(274,4 + 271,4 - 474,3) + (4 74 ,2 + 37i )2 ) a 2 ] - LE [2 ( 74 , 2 + 71,2) a + 2 73 , 3 + 4a 2 73 ,i] 

+ L 2 (2073,1 - 71,2) + 271,3 - 71,4 - [a 2 + Q] 71,2 , 
SRi = E 2 [(274,5 + 271,5 - 474,4) + (271,2 + 474,3 - 474,2 + 371,3) a 2 ] + EL [(-4 7 i, 2 - 271,3 - 274,3 + 474,2) 

x a + (473,3 - 273,4) + 473,1a 2 ] + L 2 [(— 4a73,i) + (271,2 - 71,3)] + Q (271,2 - 71,3) + 271,4 - 71,5 - a 2 7i,3 , 
6R 2 = E 2 [(-474,5 + 274,6 + 2 7l , 6 ) + (4 74 ,4 - 474,3 + 371,4 + 271,3) a 2 + a 4 ( 7 i, 2 + 274,2)] 

+ EL [(-274,2 - 271,2) a 3 + (473,4 - 273,5) + (-274,4 + 474,3 - 471,3 - 271,4) a - 4a 2 7 3 , 3 - 2a 4 73 ,i] 

+ L 2 [2073,3 + (-71,4 + 271,3) + a 2 7i, 2 + 2a 3 7 3 ,i] + Q [(-71,4 + 271,3) - a 2 7i, 2 ] + 271,5 - 71, 6 - a 2 7i >4 , (67) 

5B = -2 l3A aE 2 sm 2 6 + 2EL l3A , Se n>0 = , (68) 

<S<I>2 = E (73,1a 2 + 071,2 + 074,2 + 73,3) - 2aL"f 3A , 
f>$ 3 = E [73,4 + (74,3 + 7i,3 + 271,2) a] , 

<S$ 4 = E [(271,3 + 471,2 + 71,4 + 74,4) a + M 3 a 2 73, 3 + 73,5] - L (a 2 7i, 2 + 2073,3) , 

<5$5 = E [(71,5 + 871,2 + 271,4 + 74,5 + 471,3) a + 73,4a 2 + 73,6 - 2a 3 7i, 2 ] - L [(71,3 + 271,2) a 2 + 2a7 3 , 4 ] , 
e>$6 = E [(74,6 + 1671,2 + 871,3 + 471,4 + 71,6 + 271,5) a + 73,5a 2 + 73,7 - (271,3 + 871,2) a 3 ] 

+ L [(-71,4 - 271,3 - 471,2) M 4 - 2a 73 ,5 + 0*71,2] ■ (69) 



In the CS limit, these separated equations become 

(70) 
(71) 



5 , /MY -- 
SR CS = -M 4 f — J ELa(, 

6$ cs = -%M ( — } Ea(, 



(72) 



to leading order and 5Q CS = 0, all of which agrees with 
the results found in [2lj |. 

We finish this section with a comment on this decom- 
posed set of first order equations. The existence of an 
approximate, second-order Killing tensor allowed us to 
construct a Carter constant, which in turn allowed us to 
rewrite the second-order geodesic equations in first-order 
form. This does not necessarily imply that the resulting 
first-order equations will be decoupled (i.e., that the r 
source term depends on r only and the 6 source term de- 
pends on 8 only). Nonetheless, in the far-field limit, we 
have just found empirically in the above equations that 



the resulting first-order equations do separate. 



B. Weak-Field Expansion 

The equations presented above are not strictly weak- 
field expansions, as the GR constants of the motion 
(E, L, Q) also depend on the radius of the orbit. We 
choose to parameterize the orbit in terms of a semi-latus 
rectum p, an eccentricity, e and an inclination angle, 9t p - 
The weak- field limit corresponds to taking p ^> 1. We 
define these constants of the motion from the radial and 
azimuthal turning points of the motion, which are given 
by r(r±) = and 0(0 tp ) = 0, with 



r± =Mp/(lTe) 



(73) 



For Kerr, setting R K (r±) = = 9 K (0t P ), we find that 
^-l + ^( e 2 -l) + ^(l-e 2 ) 2 + 0(^ 3 ), (74) 
L K - v^sinfltp + 2?>~ 1/2 ( e2 + 3 ) sin6 t P 
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- (e +3) sin tp - ^ 



(3 + e ) sin tp cos 



7 t p 



sin 



tp 



8p 3 / 2 



( e2 + 3) 2 + ^(: 



2p ; 



(5 + 3e 2 ) (e 2 + 3) sin 2 tp + 0(p~ 3 ) , (75) 



2a 



p+(e 2 + 3) -- I7I (e 2 + 3) sin0 tp 



A a ■ 

4 : ,„ sm 



(e 2 + 3) + a 2 (3 + e 2 ) sin 2 tp 

tp (2 + e 2 )(3 + e 2 )) +G(p- 2 ). (76) 



p 3/2 

We distinguish here between (E, L, Q) and (-E K , L K , Q K ), 
since we are defining E K etc. to be the constants of 
the motion for the Kerr orbit that has the same turn- 
ing points as the geodesic in the deformed spacetime. 

We can now insert these relations into our expressions 
for (ST, SR, SQ, <5$), but the resulting equations are quite 
horrendous. We will thus present results only for special 
cases in which certain 7^ are no n- vanishing. If only B2 = 
(71,2,73,1,73,3,74,2) is non-vanishing (which corresponds 
to keeping only the 0(l/r 2 ) terms in h^ v ), then 

6T B2 ~ (71,2 + 274,2) M 2 - 2Af 2 a 73 ,i sin 2 , (77) 
5R B2 ~ -2A/ 2 ry/ 2 73,isin0 tp 

+ MV (71,2+ 274,2) , (78) 
5@ B2 ~ 2Af 2 p 1/2 7 3il sin0 tp - 2Af 2 a7 3 ,i sin 2 9 , (79) 
2M 3 ap 1 / 2 



(S$B 2 ~ -2 73,1 sin tp 

M 3 

+ — [a (71,2 + 74,2) + 73,3 + 73,ia 2 



(80) 



In the higher order cases, we can obtain a more general 
formula: if Bn>2 — (ji,N, 73,iV+i> 74,Jv) is non-vanishing 
(which corresponds to keeping only the 0(l/r ) terms 
in hp V ), then (for N > 2) 



ST, 



Pa- 



M 



N 



I — I (7i,jv + 274,at) 



SRb, 



rM 3 



0, 



M 



M 



N-3 



(71, N + 274, JV) , 



N 



— J [a (71, n + 74,w ) + 73,JV+iJ 



(81) 

(82) 
(83) 
(84) 



In the CS limit and to leading order in the weak-field, we 



find 



5 My/2 . 

oTcs ~ — aCsmt'tp , 

8 r 4 

5 m y / 2 . 

di? cs ~ - — aCsm^tp , 



S1> c 



0. 



5Af5 
j-aC ■ 

« 1-4 s 



(85) 

(86) 
(87) 



C. Orbital Parameterization 

If we define the orbit through the semi-latus rectum, 
eccentricity and inclination parameter introduced above, 
then in the parametrically deformed spacetimes studied 
in this paper, the turning points are not in the same spa- 
tial location as the GR geodesic with the same (E, L, Q). 
We can use the same parametrisation of the orbit if wc 
force the turning points to be in these locations. The 
corresponding constants of the motion are shifted from 
those in the Kerr spacetime with the given turning points 
by an amount of 0(e): 



E = E K + eSE , 
L = L K + eSL , 
Q = <3k + eSQ . 



(89) 
(90) 
(91) 



This only changes the (T K , R K , K , $ K ) potentials in the 
geodesic equations [Eqs. (|4T)) - (|50|) ]. as (5E,5L,SQ) cor- 
rections to (ST, 5R, SQ, 5$) will be of 0(e 2 ). 

The shifts in the energy, angular momentum and 
Carter constant can be obtained by requiring that the 
turning points of the perturbed spacetime be the same 
as those of the GR Kerr spacetime. Expanding to O(e), 
the radial turning points lead to the conditions 



= 



dR K 



dE K 



SE 



dR K 



dL K 



SL + 



dR K 



dQ K 



SQ 



SR (r± ; E K , L K , Q K ) 



(92) 



while the location of the polar turning points provides 
the condition 



= 



dE K 



SE 



dL K 



SL 



dQ K 



SQ 



+ Se(d tp ;E K ,L K ,Q K ). (93) 
These equations can be solved for (SE, SL, SQ) to obtain 



J 



SE 



R^se(e tp ) 



sin 2 6 



SR(r + ) 



tp 



^+cot 2 VR- 



R Q SO(0t P ) 
sin 2 tp 



SR(r_ 



r+ + cot 2 e tp R+) s- 1 



(94) 
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SL = 



R Q 5G{6 tp ) 



sin 



SR(r_) 



R+ - 2a 2 E cos 2 1 



't P 



sin 2 6» tp 



SQ = cot 2 6 tp 2L5L - 2a 2 E cos 2 tp 5.E - 59(0 tp ) , 
where 



V 



R, 



-4a 



(lie) 2 
P 



V 



lie 



E-2L- 



(lie) 2 
P ( V 



2a 



2 V 



E - 4a- 



lie V lie 



lie . 

Rn 



lie 



-L 



P 



SR(r + ) 



R E - 2a 2 E cos 2 6» tp i?Q J 5" 



P 



S = 



R+ - 2a 2 E cos 2 9 tp R+ 



R L + cot 2 OtpRg 



(lie) 2 lie 
R E - 2a 2 E cos 2 tp i?Q 



cot 2 6» tp i?^) , 



r 



(95) 
(96) 



(97) 

(98) 
(99) 



and all quantities are to be evaluated at (E K , L K , Q K ). 

Taking the weak-field limit of the above equations, 
when only B 2 is non- vanishing, we find 



SE-, 



S-2 



SLr 



1 

1 

VP 



(1-e 2 ) 2 (71,2 i 2a 73 ,i sin 2 6» tp ) 



(100) 



(-a sin 3 6» tp 73,i +sin0 tp (74,2 +71,2/2)) , 



(101) 

5Q B2 « -2^P sin tp 73,i, ( 102 ) 
while if only Sat is non- vanishing for N > 2 we have 

5E Bn » ip- Ar / Ar _ 2 (e)(l - e 2 ) 2 ( 7l , N + 2 74iJV ) , (103) 

« \p m ~ N fN{e) sine tp (7i,JV + 2 7 4,jv) , (104) 
<5Qb» ~ p 2_Ar /w(e) cos 2 tp (7i,iv + 2 74,a0 , (105) 
where we have define the eccentricity function 
'(l + e) w - (l-e)^' 



Me) 



(106) 



(1 + e) 2 - (1-e) 2 
If we take the CS limit, the leading order pieces are 

^cs = JaC^ 11/2 sin0tp(l-e 2 ) 2 (l + e 2 ) , (107) 



5L CS = -aQp- A sin 2 tp (e 2 + 3) (3e 2 + l) 



(108) 



5Qcs = \<P" V2 sin0 t p cos" tp {e z + 3) (3e z + l) 



(109) 



and we see that 2QSLcs = ESQ 



cs- 



IV. FUNDAMENTAL FREQUENCIES 

In this section, we calculate the fundamental frequen- 
cies of orbital motion for the parametrically deformed 
spacetime. We first compute generic expressions for these 
frequencies, and then present results valid in the weak- 
field limit. 



General Results 



1. Radial Frequencies 



The condition that the turning points of the radial mo- 
tion are at r± allows us to write the change to the right- 
hand side of the radial geodesic equation as 



<9i? Kr dR Kr dR Kc r (M\ n 

K -SE + ^SL + ^SQ + M £ SR n (- J 



8E K 8L K 
Mp 



1 - e 



3Q K 

Mp 
lie 



M 2 £ v n 



(110) 

where we have here factored out the two turning points 
r± . Equating coefficients of r™ allows us to derive a re- 
cursion relation for the constants v n : 



—V-l 



2p 



-vq 



2p 



(1-e 2 ) 

„2 



v_2 = -2E K SE , 



V-2 = , 



(1 - e 2 )" -1 (1 - e 2 ) 
2a 2 E K SE - 2L K SL - SQ + SR_ 2 



P 



-V-2 



-Vl 



2p 



(1 - e 2 ) V ° (1 - e 2 ) 
4a(aE K - L K )6E + 4(L K - aE K )SL 

+25Q + SR-! 



P 



-v-i 



(111) 
(112) 

(113) 
(114) 



-v 2 



-V n+2 + 



2p 



p 



1-e 2 ) 1 (1-e 2 ) 



2p 



: w i — ~rr — "^o = —a SQ + SRo , 

(115) 

p 2 

Vn+l ~ T, 9^ V n = &Rn Vn > 1 . 

(1 - e 2 ) 

(116) 

These equations can be solved to derive successive co- 
efficients. In practice, we will be interested in pertur- 



(1-e 2 ) 
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bations for which the series terminates at a finite N. If 
5R n >N = 0, then w„>jv = and the above equations 
provide N + 2 equations for N unknowns. The two extra 
equations are actually redundant, as they give the SE 
and SL needed to keep the turning points fixed at r±. 
These (SE, SL) were already computed in the previous 
section. 

The solution to these equations can be best studied by 
separating out a few special cases. Let us consider the 
situation where the series contains only one term, SRjy. 



Then, the above equations indicate that Vpf 



Vn<N 



P 



-2+(n-N) 



P 



and 



Hence, all terms proportional to 
v n /r n have the same leading-order behavior in p. Now, 
let us consider the following special cases: 

,v 2 ) are 

(117) 
(118) 

2 - SR-2 
(119) 



• If N < 2, then the five coefficients (v^ 2 , 
non-zero, but v n>2 = 0. 

• If N < 1, the coefficients are given by 

V-2 = -2E K 5E, 



V-l 



t'o 



4pE K 5E 
1-e 2 
P 2 (3 + 



(1-e 
P 3 (l 



2A2 



E K SE + 2LSL 



Apa 2 



(1-e 2 ) 3 
Aa(aE, 



2(2L K SL + 8Q) 
2p 



1-e 2 
L K )SE + AaE^SL 
P 



E K SE 



''2 



1-e 2 
5R-2 + SR-i 



- 1 



(5 + 10e 2 + e 4 )p 4 (3 + e 2 )p 2 a 2 



(1 



=2\4 



(1 



o2\2 



(120) 
E K SE 



(1- 
(3 + e 2 )p 2 



[AaE K SL - Aa(aE K - L K )5E] 

(2L K 5L + SQ) + a 2 5Q 



1-e 2 

p[(3 + e 2 )p - 4(1 - e 2 )] 



(1-e 2 ) 2 



6R- 



■2p 



1 



■8R- 



SR 



(121) 



where we have assumed all but one of the SR^s are 
zero. 

• If N > 1 , the previous five equations still hold, with 
the SR^s set to zero. 

• If N = 1, 2, we have the additional equations 



2p 



(1 



-v 2 



(1 - e< 

P 2 
(1 - e : 



-v 2 



SR 1: 
SR 2 . 



(122) 
(123) 



but these will be automatically satisfied for the SE, 
SL and SQ derived earlier. 



• If N > 2, we have v n — for n > N, v_ 2 , ■ ■ ■ ,v 2 
are given by Eqs. (pT7) - (fl"2T|) and 



VN-k = -(1 - e 



SR 



N 



[(l + e) fc+1 -(l-e) 



fe+11 



2ep 2 + k 

for < k < N - 2 . 



We note that Vk is proportional to p k ( 2+N \ 

The set of equations for the w n 's, SE and SL are linear, 
so to find the general solution we can proceed as follows. 
We denote by v~ 2 , SE_ 2 and SL_ 2 the solution to the 
above equations with only 5R- 2 ^ and define 

SQ„ 2 = cot 2 8 tp 2L K 5L- 2 -2a 2 E K cos 2 tp 6E- 2 -6Q(9 t p). 

(124) 

Then, for N > —2, wc denote by v% , SEn and SLn the 
solution to the above equations with only S R n ^ and 
with 

SQ N = cot 2 6 tp 2L K SL N - 2a 2 £; K cos 2 8 tp SE N . (125) 

The general solution is then given by v n = ^2 k v k , SE = 
J2 k SE kl SL = J2k SL k and SQ = J2 k 5 Qk- 

The radial geodesic equation is most conveniently in- 
tegrated by changing variables. First, we define a new 
dimcnsionless time parameter, A: 



dX 



(126) 



where r is proper time (differentiation with respect to 
which was denoted by an overhead dot previously). Next, 
we parameterize the orbit in terms of the (dimcnsionless) 
semi-latus rectum p, the eccentricity e and a phase angle 
if) via 



Mp 



1 + e cos ip 



(127) 



Equation (jllOp then becomes a differential equation for 
ip, namely dip/d\ = y / V^, where = Vki/> + cSV^ and 



V K ^(ip,P,e,i) 



SV^(ip,p,e,i) 



1-Ei 



(1-e 2 ) 3 

w 



[p(l - e 2 ) -p 3 (l + ecos^)] 

-p 4 (l+eCOS1p)] , (128) 



n+2 



e 2 ) ' — ' p 

1 n=-2 y 



(129) 



The Kerr potential depends on Mp 3 = (1 — e 2 )r 3 and 
Mp4 = (1 — e 2 )r4, where r 3i 4 are the other two turning 
points of the radial motion. These are given in terms of 
E K , Q K , p and e by the expressions 



P3 = Y+ ^Y 2 - X, 
where 

1-e 2 



Pi 



= Y — \JY 2 -X (130) 



Y 



1-E 2 



P, 



X 



(1-e 2 ) 3 a 2 Q 
1 - E 2 v 2 



(131) 
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These relations imply that p 3 ~ 1 and P4 ~ a 2 in the 
weak field, so V^, « (1 — i? 2 )p 2 /(l — e 2 ) when p>l. 

We can integrate the evolution equation for ip over a 
complete orbital cycle to determine A r , i.e., the A time 
elapsed per radial cycle. Writing this in terms of the 
difference to the Kerr value, A r = Ak t + eSA r , we obtain 
the correction to the radial period 



SA r = 



2tv 



dtp 



2ji 



\/Vi> Jo 
(1-e 2 ) 7 / 2 



(1-^)3/2 



„P Jo 



ecostp) 



n+2 



x [p(l — e 2 ) - p 3 (l + e cos ip) 



-3/2 



x [p(l-e 2 )-p 4 (l + ecosV')] 3/2 d^ , (132) 
which is valid for arbitrary radius at leading order in e. 

2. Polar Frequencies 

If we write z = cos 2 6, the condition that the turning 
points of the polar motion are at z = cos 2 9 tp allows us 
to write the change to the right- hand side of the polar 
geodesic equation as 



(1-z) 



[SQ - z{5Q + 2L K SL) + 2a 2 E K SEz(l - z)] 



2M 2 73 ,! [E K L K - aEl{\ - z) 



M 



, cos 2 9 



t P 



(to + hz) 



1 - z 

and by matching powers of z we find 
1 



(133) 



in 



cos^ 



7tp 



[5Q + 2 l3 . 1 {E K L K ^aEl)} , (134) 



h = 2a(aE K SE + l3A El). 



(135) 



The coefficient of z then gives the relationship between 
8Q, 8E and SL we found in Eqs. (|94 |) - ([96 ]) . 

Let us now rewrite the polar geodesic equation in a 
simpler way, by parameterizing the polar motion using 



cos 2 9 = cos 2 9 tp cos x ■ 



(136) 



The polar equation then becomes 



d% 
dA 



Vk x (x,P, e 7 + e(*o + h cos 2 6» tp cos 2 \) , 

am 



VK x (x,P,e,t>) = a 2 (l - El)(z+ - cos 2 9 tp cos 2 x) , 

Qk 



where z + = 



(1 - £ 2 )cos 2 i 



'tp 



(138) 
(139) 



As in the radial case, we can compute the perturbation 
to the 9 period in A by writing Ag = Ak g + eSAg and 
find 

1 r h+ h cos 2 6» tp cos 2 x 
SAg = — / —dx ■ 



(l-E 2 ) 3 / 2 J Q ( z+ -cos 2 tp cos 2 x) 3/2 



(140) 



3. Azimuthal Frequencies 

We can write the azimuthal geodesic equation as 
d<p p 2 d(p 



dA M dr 



t> r (r)+<5> e (9), 



(141) 



which allows us to define a radial and a polar contribution 
to the average advance of (p over a radial /polar period 

A <^r = / 1 ' \ . 7TT —d^ , (142) 



o d^/dA 
27r ^[cos-^cosgtpcosx)] 



II 



dx/dX 



dx- (143) 



The average rate of advance of <f>, w^a, can then be com- 
puted as 

W M = ( — ) - — — + — — ■ (144) 



dA/ A A r 



A e 



We can write the individual contributions in the form 
A(f>A, r = AcpKA.r + eSA(f>A,r and A(pA,e = A(p K A,e + 
e5Acf>A.e in which we take the Kerr pieces to be 



AcpKh.e = —TP, 
sin 



(145) 
(146) 



Each of the e corrections has two contributions — one 
from the change in the numerator, $ r etc., and one from 
the change in the denominator, dtp/d\. We obtain 



<5A0A,r = a 



(1-e 2 ) 3 / 2 f 27r 2p{l + ecosip)5E- a(l + ecosi>) 2 5L 
(1 - El) 1 / 2 J D r {i))[p 2 - 2p(l + ecostp) +a 2 (l + ecostp) 2 ] 



dtp 



14 



(1 - e 2 ) 3 / 2 ^ s^ n f 27T (l + ecosVQ" ; 
(!-^) 1/2 ^ P n Jo 



DM) 



(1-e 2 ) 7 / 



2 00 



ff (1 + ecos-0)™ +2 {2p(l + ecosip)E K - a(l + ecosip) 2 L K } 



(1 - £ 2 ) 3 / 2 ^ p n J DMT [p 2 - 2p(l + ecos^) + a 2 (l + ecos^) 2 ] 

DM) = { [p(l - e 2 ) -p 3 (l + ecos^.)] [p(l - e 2 ) -p 4 (l + e cos </>)]} 1/2 , 



(147) 
(148) 



and 



£A</> A ,e = 



1 



2;r 



a-v/ 1 Jo (1 - cos 2 6> tp cos 2 x) a/ z + - cos 2 tp cos 2 x 
1 fa + ji cos 2 gtg cos 2 x)-^k 



dx 



a 3 (l E 2 ) 3 / 2 J (1 — cos 2 6>t p cos 2 x) (^+ — cos 2 6* tp cos 2 x) 3 ^ 2 



dx- 



(149) 



In both the radial and polar corrections, the first term 
is the contribution from the change in the orbital con- 
stants E and L in the Kerr part of the numerator. The 
second term is the contribution from the change in the 
numerator due to the addition of the perturbation to the 
numerator. The third term is the contribution from the 
change in the denominator. There is no second term in 
the 9 contribution since the perturbation to <& is purely 
radial. 



4- Temporal Frequencies 

The motion in t can be treated in the same way as the 
motion in <f>. With the analogous definitions, the average 



I 

rate of advance of time with respect to A, ujta, is 
dt\ AT A . r AT A , e 

ujta - 



dX/ A ' A fl ' 



(150) 



where AT\^ r = ATk A,r + eSAT\ tr and AT^g = 
ATk a,0 + £5ATa,8- The Kerr pieces of this rate of change 
are 

AT K A,r = (r 2 + M 2 a 2 )^, 

AT^A.e = -aM(aME K sin 2 9 - ML K ) , (151) 
while the 0(e) corrections are 



5AT Ar (1-e 2 ) 3 / 2 f 27r [p 2 + a 2 (l + ecosV') 2 l {[p 2 + a 2 (l + ecos^) 2 l<5£;-a(l + ecosV) 2 <5L} , 

— -dip 



M (l~E 2 ) 1 / 2 J (l + ecosi/0 2 -Dr(V0b 2 -2p(l + ecosV') + a 2 (l + ecosV') 2 

(1 - e 2 ) 3 / 2 SJZ [ 2 * (1 + ecosVQ" , , 
(l-^^p-io DM) V 

(l_ e 2)7/2 ^ v„ f v (l + ecosV) n [p 2 + a 2 (l + e cos t/>) 2 ] 



,_ J JO 



1 - £ 2 ) 3 / 2 p n J DM) 3 [P 2 - 2p(l + ecos?A) + a 2 (l + ecosV>) 2 ] 



and 



x {E K [p 2 +a 2 (l + ecosi/') 2 ] - aL K (l + ecosip) 2 } dip , (152) 



SAT Ae 1 f 2 " a5L-a 2 (l - cos 2 9 tp cos 2 X )SE , 2 73 X E K f 2 " (1 - cos 2 9 t „ cos 2 x) , 

- d X 7==5F / / F= o d X 



M a-v/1 - 7o vi+ - cos 2 6> tp cos 2 x y / T^E 2 Jo \/z + - cos 2 9 tp cos 2 x 

1 /" r (*o + *i cos 2 6> tp cos 2 x) [aL K - a 2 E K (l - cos 2 9 tp cos 2 \)] 



a 3(l_ ^2)3/2 y Q (z + -COS 2 tp COS 2 x) 3/2 ' (153) 

I 

in which D r (ip) was defined in Eq. (|149[) . The origin note that 5T£ denotes the radial part of the perturbation 
of each term is as described in the azimuthal case. We to the temporal potential. This is equal to 5T n , except 
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STS = (71,2 + 2 74 , 2 )M 2 



5. Combining Results 

The orbital frequencies quoted so far were written in 
terms of the A parameter. We compute the perturbations 
to the physical frequencies, i.e., those expressed in Boycr- 
Lindquist time, by first defining 



AT A = A AT A , r + A r Ar Alfl 

A0 A = AgA(j)\ r + A r A0 A a 



(154) 
(155) 



and then writing 



iljy — 



n r = 



2-nAg 
AJ\ 



2ttA,. 



(*) AT A ' AT\ AT A ■ 

(156) 

The factor of 2ir is included in the definitions of fl r and 
Qg so that all three frequencies are angular frequencies, 
expressed as radians per second. From the preceding 
equations, we may find the perturbations to the physical 
frequencies to be 



SCI4 



5AT A 



A<j) A AT A ' 

SOg _ SA r 
fig A r 



SAT A 

TT = aTF". (157) 



5fl r 5A 

n r — a 

SAT A 



AT, 



A 



(158) 



The frequencies that we need for the construction of 
AK waveforms, using the notation of Barack and Cut- 
ler (l^ | . are the orbital frequency, 2ttv, the perihelion 
precession frequency, 7, and the orbital plane precession 
frequency, a (see Section [VI[) . These are given by 

2ttv = £l r , 7 = Cl$ — O r , a = f^0 — 0,0 , (159) 

and therefore the changes in the frequencies can be com- 
puted from 



5v _ SA e 5AT A 

~~Tg~ ~ AT A ' 
S^ 6{A r -A e ) SAT A 
if Ar-Ag AT A ' 

a _ 5(Act) A - 2nA r ) SAT A 



a A(f> A — 2-kA t 



AT a 



(160) 
(161) 
(162) 



B. Weak-Field Expansion 

1. Radial Frequencies 

Denoting by v& N i the leading order in p piece of the 
parameter Vi under assumptions Bn, we find for B2 

(1 - e 2 ) 2 

V B 2 , -2 = ^3 [71,2 + 2asin 2 6» tp 73,i] , (163) 



(1 _ e 2) 

^82,-1 = -2 V — ; [71,2 + 2a sin 2 tp73il ] , (164) 



v B 2 ,a 



p* 
P 



[71,2 - 2a sin 2 0tp73,i] 



(165) 



and the other vs 2 i 's are subdominant in the sense that 
v B 2 ,i/P l 1S higher order in 1/p. We see that for these 
dominant terms vb 2 t/p 1 ~ 0(1 /p), which motivates us 
to define p- independent terms vb 2 f = v& 2 i jp l ~~ x . For Bn 
with A?" > 2, we find v& N Jp l ~ 0(l/p w ~ 2 ) so we can de- 
fine the p-independent quantities Vb n ( = Vb n i /p' l+2 ~ N . 
We find vb n { ~ for i > N — 4, 

%„,_ 2 = -(1 - e 2 ) 2 / A r_ 2 (e)(7 1 , Ar + 274,jv) , (166) 

and for -1 < i < N - 4 

v BNii = -2(1 - e 2 )/ A r_ l _ 3 (e)(7i,JV + 2 74 ,w) , (167) 



in which /jv(e) is as defined in Eq. (|106[) . 

The corresponding expressions in the CS limit are 



^CS,-2 


= --(1- 
2 V 


e 2 ) 2 (l + e 2 )Casin0 tp , 


(168) 


"CS,-1 


= -5(1- 


e 2 )(l + e 2 )Casin6» tp , 


(169) 




= --(1- 

4 V 


e 2 )(3 + e 2 )Ca sin6> tp , 


(170) 


"CS,1 


= --(1- 
2 V 


e 2 )Ca sin 9 tp , 


(171) 


"CS,2 


= --(1- 

4 V 


e 2 )Ca sin 8 tp , 


(172) 


in this 


case v C s N 


„=vcs N ,Jp l - 7/2 - 





With this identification of the leading-order parts of 
the Vi's, the integral for the change to the radial fre- 
quency reduces to 

00 

SA r = -P~ K (1 - e 2 )- 1 J2 V n In+2 , (173) 
n=— 2 

in which J„ = f (1 + e cos^)"^ , (174) 
Jo 

the v n 's are as defined above and K represents the scaling 
of the frequency correction with p, which is K = 5/2 for 
B 2 , K = N - 1/2 for B N with AT > 3 and K = 5 in 
the CS limit. We note that although we have written oo 
as the upper limit of the sum for convenience, in general 
we will be able to terminate the summation as described 
above, e.g., at n = for £>2, at n = N — 4 for Bn with 
N > 3 and at n = 2 for CS. The integral I n is evaluated 
in Appendix [X] and shown to be equal to 



,l= ^^l2/ {n-2k)\{k\f 



(175) 



in which [^J denotes the largest integer smaller than x. 
For the various special cases we have been considering 



16 



elsewhere, the leading-order shifts in the radial frequency 
can be simplified to 



SA^ 



22,5/2 



[(4-3e 2 ) 7l , 



2a (8-e 2 ) sin 2 t p 73) i] 



6A^ = 



2p 5/2 



B 4 _n(8-e 2 ) 



5 A** = 



5 A? 5 = 



2pV 2 
tt(15 + 5e 2 - e 4 ^ 



(7i,3 + 274,3) , 



(71,4 + 274,4) , 



(71,5 + 274,5) 



(176) 
(177) 

(178) 

(179) 



8A c r 



5n 
3V 



a< sin tp (12 - e 2 ) (8 + 9e 2 ) . (180) 



2. Polar Frequencies 

For p > 1, we find z+ - (Q + L 2 )/a 2 (l - £ 2 ) and 
Q + L 2 w p, so Vk x (x,p, e, t) p. The weak-field limit 
is therefore 



P 



,3/2 



(181) 



which for assumptions B2 gives 

-jj^ {71,2 + 274,2 + [3acos 2 6> tp - 4a] 73,1} , 

(182) 

for assumptions Bn with N > 3 gives 

<$A e£!jv w - N _ 1/2 fN{e) (7i,jv + 274,at) , (183) 
and in the CS limit gives 

7T 
4 p 5 



-^aCsin0 tp (3 + e 2 )(l + 3e 2 ). (184) 



3. Azimuthal Frequencies 

In the weak-field, the leading order part of <5A0A, r 
comes from the terms 



<5A0 A ,r 



AitaSE 



OO r- -,- L 2 J 

(5$, /PA 2fc 



p 



3/2 



-vfE^E(l) 

7 — AT r 1 n 



l=N^ 1 k=0 



II 



2a 



(Z-2fc)!(fc!) 2 (l-e 2 )p K + 



J V n In+3 



(185) 



in which we are using to denote the first non-zero 8<f>i, 
I„ is as defined above and K has the same meaning as in 
the radial case: K = 5/2 for B 2 , K = N— 1/2 for B^ with 
N > 3 and A' = 5 in the CS limit. Under assumptions 



B2, the second term can be seen to dominate, N, 
and 8<&2 ~ — 2a-»/p sin #tp 73,1 • We therefore find 



<5A<^ r = -^073,1 (e 2 + 2; 



(186) 



Under assumptions Bn with N > 3, the second and third 
terms both contribute and we find 



<5A^™ >2 = N ]_ 1/2 {^[l3,N+i + a(7i, w + 74,jv)]ijv 

JV-l 

+ 2a 7r(l-e 2 )/Ar_ 2 (e)+2^ f N -k(e)I, 

k=2 

x ( 7i ,jv + 27 4 ,at) } , (187) 
which for the first few TV's gives 



1 



5A ^r = ^TJ (( 2 + 3e 2 )73,4 + o(5 + 3e 2 ) 7 i, 3 
+ 0(8 + 3e 2 ) 74 , 3 ) , 

5A <:^^(( 8+24e2 + 3e4 )73,5 

+ a(40 + 36e 2 + 3e 4 )7i, 4 
+ 3a(24+16e 2 + e 4 )74, 4 ) , 

5A4% r = ((8 + 40e 2 + 15e 4 ) 73 ,6 

+ 2a(34 + 50e 2 + 9e 4 )7 1 , 5 
+ a(128 + 160e 2 + 21e 4 )7 4 , 5 ) ■ 



(188) 



(189) 



(190) 



In the CS limit we find that the second term again dom- 
inates, which gives 



• r )7T 



a( 1 + 3e z 



The dominant contribution to 5A(f>\g is given by 



(191) 



<5A0 



9 3/2 



2Q K 8L - cos 2 tp t o L K 



cos^ W tD sin 



t p sm utp 



SL + ^Lk 1 



2a 2 (l-e 2 ) 



With the #2 assumptions, we find 

SA(j)® 2 e = — sin6» tp 73,i 



7 t P 



(192) 



(193) 



For assumptions Sat with TV > 3 and in the CS limit, we 
have z_irj = <5Q and using Eqs. ([94 |) - (j96]) , we see that 



2Q K SL - L K SQ 



s 2 9 



t P 



Qk 



i 2 - 



s 2 6 



t P sin tp 



2a E K L K 5E , 
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2a 2 
P 



(l-e 2 )SL + 2a 2 E K L K SE. 



(194) 



Both terms therefore contribute at the same order giving 



h B N>2 _ Q 2 tt(1 - e 2 ) 

9 A,0 _ p N+ 



^±[f N (e) + (l-e 2 )f N - 2 (e)] 
x (fi,N + 274,w) sin tp , (195) 



which for the first few iV's gives 



<5A0 



2a 2 7r(l - e 2 ) 



3 (7i,3 + 2 7 4, 3 ) sin tp , (196) 



5A<f£ a 



3a 2 ^(l-e 2 ) / e 2 \ , „ , . „ 

*e = - 5 " ( 1 + 7 ) (^,4 + 2 74,4) sm0 tp , 

(197) 
? • 

(198) 



_ . h 4a 2 7r(l — e 2 ) , , . 
5A^ e = 4 i(l + e 2 ) ( 7li5 + 274,5) sin0 tp 



In the CS case, we find, at leading order in p, 



5A(/> 



25tt a 3 C 2 , 



A.e 



4 p 13 / 2 



1-e/ 1 + 2^ + — sin 



(199) 

If we linearize in a, for consistency with the order to 
which the CS limit was derived, this term can be seen to 
vanish SA^g = 0. 



4- Temporal Frequencies 

In the weak-field, the leading order part of SAT\,, 
comes from the terms 



SAT A , r 2^I 2 SE ~ 8T1 ^ 



M 



(l_ e 2)3/2 ■ -v- ^ (+1 ^ V2 

Z=JV T 



/! 



fc=0 
AT-4 



2A- 



(Z-2fc)!(fc!) 2 (l-e 2 )p 



2\ rt K-2 



(200) 



n=-2 



in which Nt denotes the first non-zero 5T n and I n , K 
have the same meaning as in the radial case and as in 
the preceding section. 

We can now evaluate the dominant contribution to 
<5ATa j7 . in the weak-field for the various cases we have 
been considering. For assumptions B2, the second and 
third terms dominate and wc find 



Under assumptions Bn for N > 3, the third term domi- 
nates. For the first few iV's we obtain 



SAT t = 2p l/ 2( 3 1 7r _ e 2)l/2 ( 2 T4,3 + 71,3) 



(5 AT: 



ttM 



P 



3/2 



SATt = *M 



(l_ e 2)i/a 
3 + e 2 

(l_ e 2)l/2 



(202) 

(2 74 ,4 + 71,4) , (203) 
(2 7 4,5 + 71,5) ■ (204) 



In the CS limit, at leading order in p, the third term 
again dominates and we obtain 



157dl/Ca . „ 

— 3 sin0 t p 

8 p A 



4 (! + e2 ) ,1 2^ 



(205) 



The dominant contribution to <5ATa,0 comes from the 
terms 

SAT A: e _ 2ira8L 2ira~/ 3 . 1 E K (2 - cos 2 tp ) 



M 



VP 
TraL K 
p3/2 

7ra 2 T K 



Vp 



to 



cos tp ti 



P 



3/2 



IS 2 



tp 



1 - 



3 cos 2 tp 



to 



cos 9 tp ti 



(206) 



Under assumptions B2, the second term dominates and 
wc find 



SAT 



AM 



M 



2ira 

Vp 7 



73,1(2- cos 2 6 tp ). (207) 



Under assumptions Bm for N > 3, the final term domi- 
nates giving 

5AT® N e >2 7r a 2 

" = ^73172-/^ ( e ) (2 -cos tp ) (274, N + 7i,jv) ■ 



M 



(208) 

At leading order in p the same term dominates in the CS 
limit to give 

6AT$ 5tt a 3 C 



5 (2-cos 2 0tp)sin0 t p (e 2 + 3) (3e 2 + l) . 

(209) 



If we linearise in a for consistency, then wc find 5AT£ s g 
in the CS limit. 



5. Physical frequencies 



SAT®*. = ^ [ 7l , 2 + 4 74 ,2 - 2a sin 2 tp73l i 



(1-e 2 ) 1 / 



1/2(71,2 + 2a sin 0t P 73,i) 



(201) 



In the weak-field Kerr metric, the leading-order parts 
of the frequencies are given by 



2tt 2tt 

A r if — —- JA0K — —=L , 

Vp Vp 



(210) 
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AT A ^ rK 

M 
M 



2na 

—— A4>A 0K = 27T . 

VP 



= 2tt 



P 



1 -e 2 
27ra sin tp , 



l-e' 



AT AK 
A rK — A SK 



4n 2 Mp(l - e 2 )" 3 / 2 ; 
6irp- 3/2 , 



A(f)\K — 27rA r j^ = 87r 2 ap 2 
from which wc obtain 
3(1 -e 2 ) 3 / 2 



(211) 

(212) 
(213) 



(214) 
(215) 
(216) 
(217) 



,5/2 



(1 - e 2 ) 3 / 2 
Ma K = 2a± ^ . (218) 



We can now put together the pieces from the preceding 
sections to derive the leading order corrections to the 
three frequencies, which wc require in order to construct 
the modified AK waveforms. 

a. Orbital frequency. The perturbation to the or- 
bital frequency, 2-kv, is 

M6(2wv) = MSQ r = MCl r - 

\ Ag AT A / 

(l-e 2 ) 3 / SATa.t , oAT A , 



2np 5 / 2 \ My/p Mjp 
+ asin(9 t p [5A r - 5A 8 ]) . (219) 

Under assumptions Bi , the first two terms dominate and 
wc find 



(l-e 2 ) 3 
2p 7 /2 



1 



7i,2 + 47 4/ 



2a 1 



4- 



sin ftp 73,i 



(220) 



Under assumptions Bjy with N > 3, the first term alone 
dominates and we find 

MSn " N>2 = hp-^TT 2 E w - (221) 



n=-2 



For the first few N's we have 
3(1 -e 2 fl 2 



Bs _ 



4 p 7/2 

(1 - e 2 f' 2 
2p s / 2 

3(1 -e 2 ) 5 / 2 
4pii/2 



(2 7 4,3 + 7i,s) , (222) 
3 + (l-e 2 ) 1 / 2 ! (274,4 + 71,4) , 



(223) 



3 + e 2 + 2(l-e 2 ) 1 / 2 



x (274,5+71,5) 



(224) 



In the CS limit, the first term again dominates and, lin- 
earizing in a for consistency, we have 



so 



CS 



15M(l-e 2 ) 5 / 2 



16p e <a sin dt P 



4(l + e 2 ) + (4 + e 2 )(l-e 2 ) 1 / 2 . (225) 



b. Perihelion precession frequency. The perturba- 
tion to the perihelion precession rate, 7, is given by 
Eqs. f|160[) - (|162[) . which can be simplified to 



= (A 9 SA r - A r 6A 6 )(AT Ar + AT Ae ) 
7 (A, - A g )AT A 

_ (A r SAT Ae + A e SAT Ar ) 
ATa 



(226) 



In all the cases we have been considering, the first term 
makes the dominant contribution in the weak-field and 
wc find 



(l-e 2 ) 3 / 2 
MSj = 1 _ ; (SA r - SA e ) . 

ZTTp 



(227) 



Under assumptions B 2l only the SAg term contributes 
giving 



M _ p2\3/2 



2p 5/2 [T1.2 + 274,2 
3a cos 2 6» tp - 4a) 73,1] . 



(228) 



Under assumptions Bn with N > 3, both terms con- 
tribute and we find 

MS ^n>* - (1 ~ e2)3/2 



2pN+ V2- [f N (e) + {l-e 2 )f N . 2 {e) 



N-2 



+ 2j2fN-k-i(e)h 



k=l 



(71, AT + 274,at) , 



which for the first few N's gives 



B ,_ 3(l-e 2 ) 3 / 2 
2 P V 2 



(71,3 + 274,3) 



(229) 



(230) 



Qf1 _ p 2\3/2 

M ^ 4 = V /2 ( 4 + e X^M + 274,4) , (231) 



, Bb 5(1 -e 2 ) 3 / 2 
Ap 1 



T7 f-(4 + 3e 2 )( 7l , 5 + 274,5). (232) 



In the CS limit, both terms again contribute and we find 



M67 



75(1 -e 2 ) 3 / 2 
64p 6 



aCsin6i t p(8-|-12e 2 + e 4 ). (233) 
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c. Orbital plane precession frequency. Equa- 
tions (|160p - (|162j) for the orbital plane precession 
frequency can be written as 

8a _ A,.(5A0 e + A e 8Acj) r + 8A e A<j) r + (A0 e - 2n)SA r 
a 



A0a - 2nA r 



6AT A 
AT A ■ 



(234) 



In the Kerr limit, A(f)g — 2ir ~ 1/p 2 . Using the preced- 
ing results we can see that in all the special cases we 
have considered in this paper, the term on the first line 
dominates over the term in the second line. For B2, the 
dominant contribution comes from the 8A<j)e term 



MSa® 2 = 



/■1 _ p 2\3/2 

p5 /2 osmgtp73,i- ( 235 ) 



If 73.1 = 0, the B2 limit instead gives 



M8a L 



a{l-e 2 f/ 2 



P 



(71,2 + 274,2) • (236) 



For assumptions Bn with iV > 3, the dominant contri- 
bution comes from the SAa term and we obtain 



M5a L 



i(l - e 2 ) 3 / 2 



2p N + 



J —fN{e)(ji,N + 2j 4)N ), (237) 



which for the first few Ns gives 

MSa Bs = ~(1- e 2 f 2 (3 + e 2 ) (71,3 + 2 7 4, 3 ) , 



(238) 



M5a l 



(l-e 2 ) 3/2 (l + e 2 ) (7M + 2 74 ,4) , 



(239) 

= ~ (1 - e 2 ) 3/2 (5 + 10e 2 + e 4 ) (71.5 + 2 74 , 5 ) 



(240) 



In the CS limit, the dominant contribution comes from 
the 8A<f> r term and we find 



MSa c 



5 Cl - e 2 ) 3 / 2 
a C J (8 + 24e 2 + 3e 4 ). (241) 



G4 



d. Kepler Law in the Circular Limit. With all of 
these frequencies, it is instructive to compute the cor- 
rection to the Kepler orbital frequency when considering 
orbits in the equatorial plane. This frequency is given by 
f^Kcp = Q<t> = 27w + a + 7, and for GR this is simply 



MQ Kt 



P,GR 



P 



3/2 



(242) 



For the Kerr deformed metric, fl Kcp = l] Kl , p , GR + edfl Kcp , 
where in the Bn limits 



^Kc P ,e 2 



n5/2 



(71,2 + 274,2 + 073,1) (243) 



5CI 



K0P.S4 



Kcp,£?5 



3 7 i,3 + 274,3 

4 pV 2 
7i,4 + 274,4 

p 9/2 

5 7i,5 + 274,5 



4 p»/2 

and in the CS limit 

CS 



5 Co 
4p6 



(244) 
(245) 
(246) 

(247) 



which agrees with [4l|. Notice that this correction is of 
4.5PN absolute order, which is consistent with the fact 
that the CS correction to the Kerr solution enters at 3PN 
order relative to the 1.5PN order spin-orbit coupling in 
the gravitomagnctic sector. 



V. IMPLICIT DEFORMATION OF 
RADIATION-REACTION 

We have so far concentrated on corrections to the con- 
servative sector of GR, i.e., to the metric and the shape 
of orbits, but have not yet considered dissipative cor- 
rections, that is, modifications to the radiation-reaction 
force and the fluxes of energy, angular momentum and 
Carter constant. Modifications to the orbital shape will 
implicitly translate to modifications to the fluxes, as the 
latter are computed from derivatives of the quadrupole 
moment, with the particle moving on a geodesic orbit as 
a source. In what follows, we will compute these implicit 
modifications to the radiated fluxes. 

We will, however, neglect any direct modifications to 
the fluxes. Stein and Yunes [31( have shown that there 
is a wide class of quadratic gravity theories for which the 
effective gravitational energy flux is indeed still described 
by the Isaacson effective stress-energy tensor. An exam- 
ple of such a theory is dynamical CS gravity [l8| . There 
can be, however, theories in which dynamical scalar fields 
emit dipolar radiation, which in turn might lead to mod- 
ifications of similar order as those considered here. We 
leave the inclusion of such direct dissipative modifications 
to future work. 



A. Equations of Motion 

To estimate the implicit corrections to the radiated 
fluxes, we must evaluate the leading order in e and p 
correction to the quadrupole moment tensor. This can 
be accomplished by linearizing the geodesic equations in 
e and then separately linearising the Kerr geodesic part, 
i.e., the O(e ) part, and the linear correction, i.e., the 
0(e) part, in p. Using the same orbital parameterization 
as before, in terms of the three phase angles (ip, \, 4 > )i we 
find that in this approximation 



d^ = ^ (1 + eC0S ^ 



20 



1 



2(l-e 2 )p K + 



T ^"^ + ecos V0 



n+4 



SE 
pT/2 

STn t 



P 



(1 + ecosi/") 2 

(l + ecosV0 WT+4 



N T +7/2 



dx 
dt 



^T = 372( 1 + ecos V') 2 



t\ COS 2 #tp cos 2 x 

2p~^ 
SE 

-3/2(1 + ecos^) 2 



2p 5 / 2 
(1 + e cosf/*) 2 



(248) 
(1 + ecosVO 2 



■(1 + e cos 



JVt+4 



d0 

dt 



pAT T +3/2 

sin^tp (l + ecos?/;) 2 

p3/2 I _ cog 2 CQS 2 ^ 



(249) 



-(1 + ecos?/')' 1 



oX (l + ecos'0) 2 
p 2 1 — cos 2 8 tp cos 2 x 

^(l + ecos^ +2 
o\E , sin# tp (1 + ecosip) 2 

p3/2 l _ cog 2 cog 2 ^ 

(5Tjy T sin6> tp (1 + ecos^) jVr + 4 

pN T +7/2 I _ cog 2 ^ p cog 2 ^ 



The first term in each of these equations is the Kcplerian 
solution and (Nt, N$) denote the leading-order terms in 
the perturbations 6T and (5$ as before. The terms in- 
side the brackets are of different orders in p in general, 
but the dominant terms differ under the various sets of 
assumptions we have been considering. 

Under assumptions $2, the first and third terms are 
dominant in dip/dt, the first two terms are dominant for 
d\/dt and the second term is dominant for d<p/dt. Our 
approximation to the orbit then becomes 



dtp 1 



1 



df 



,3/2 



(1 + ecosVO 



,7/2 



71,2 + 274,2 (1 + e cos ip) 



(1 - e 2 ) 

+ (1 + e cos V>) 71,2 + - — s — -(1 + e cos VO 71,2 



2a 



1 



(1-e 2 ) 



(1 + e cos ip) 4 + (1 + e cos ip) s 



+ -(1 + ecosVO" 7 3 ,isnrw tp 
+ 2a(l + ecos'0) 4 73,i(l — cos 2 6? tp cos 2 x) 



(251) 



dv 82 1 1 

— = pW ( J + e cos ^ + 6 K 2 + 2 ^. 2 

+ 2(cos 2 6» tp - 2)073,1 + 2(173,1 cos 2 tp cos 2 x) 
x (1 + ecos^) 2 , (252) 

1 /7l,2 



icp® 2 _ sin^tp 
~dT ~~ p 3 / 2 l-cos 2 6» tp 



(1 + ecosip) 2 
? 6 tp cos 2 x 



'74,2 



(l + ecosi/0 2 , . 
asm fl^)^ —-———- . (253) 

Under assumptions Sat, the first term dominates in 
dip/dt and in dx/dt, while the second term is dominant 
in d<p/dt. Our approximation to the orbit then becomes 



dV> 



|SjV>2 



1 



df 



■(1 + ecosip) 2 



p3/2" 

/ ^^/A r -2(e)(l + ecos^ 2 



6 yV+l/2 



W-4 

/jv-i-3(e)(l + ecosi/') Z 

i=-i 

(71, AT + 274,jv) 



X 



(254) 



-^ = ^37T( 1 + ecos ^) 2 

-f N (e)(l + ecosip) 2 (7i,jv + 27 4) iv) , 

(255) 



E 2p w+i/2- 



(250) d(/) Bw > 2 _ sin6> tp 



(1 + ecosV>) 2 1 , , s 

p 3/2 I _ cos 2 Q tp cos 2 x + 6 2p N+l/2 JN ^ e > 



dt p 

(1 + ecosi/') 2 

(7l,iV + 274 N ) r-Tj 5— 

1 — cos 6>tp cos J x 



(256) 



In the CS case, the first term dominates in dip/dt and in 
dx/dt, while the second and third terms are dominant 
d(p/dt. Our approximation to the orbit then becomes 



dip' 



,cs 



5 

e cos ip) 2 — e— -a( sin 9 tv 



^T 2 

X [2(l-e 2 )(l + e 2 )(ln 
+ 4(1 + e 2 )(l + 60081/0' 
+ (3 + e 2 )(l + ecosV') 4 
+ 2(1 + ecos?/') 5 + (1 + ecosi/0 6 ] , 
d X cs 1 ,1 



dt p 

- e cosip) 



^aCsin0 tp (3 + e 2 )(l + 3e 2 )(l 



i(p cs sin0 tp (1 + ecostp) 2 
dt p 3 / 2 1 — cos 2 0t p cos 2 x 

+ e ^ a C^(3 + e 2 )(l + 3 e 2 ) 



(257) 



3 cos?/)) 2 , 
(258) 
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. o (1 + ecosVO 
x sin' tp - 



1 1 — cos 2 9 tp cos 2 x 



Fluxes 



(1 + ecos^) e 



(259) 



Using these simplified equations of motion, we can 
now compute the corrected fluxes. For this, we assume 
that [S M 



1 

5// 



E) = -—( UjT' 



1 

- I ii I jj 



Li 



2 / -j —kr, 



(260) 



where p, is the reduced mass, E and Li are the dimen- 
sionless energy and momentum defined in Eq. (|36[) . the 
angle-brackets stand for an average over several wave- 
lengths and we have defined the quadrupolc moment 



/" // •'(/; z j (t). 



(261) 



in which the particle trajectories are z % = 
[r cos <psm9, r sin</>sin#, r cos 9]. In the AK model, 
it is assumed that the orbital inclination remains con- 
stant, as the change in inclination is a higher order PN 
effect. If we make the same assumption here, 9 tp — 0, we 
can find the rate of change of Carter constant in terms 
of the rates of change of energy and angular momentum. 
This gives 



2^/p cos 2 9 tp 



sin 9 



5L - 2a 2 cos 2 9 tp SE 



„ 2 2o fdSE . dSE . 
2a cos tp [-q^-Pk + ~^r e K 

2 v /pcos 2 6> tp fdSL . 35L . 

— - n, r A p 

sin6> tp V dp 
fdSQ. dSQ . 



Pk H ^ — c K 



V dp 



■Pk 



de &K 



(262) 



In all cases except B2, at leading order this is equivalent 
to writing 



Q 



(263) 



where, as before, an L without a suffix is being used 
to denote the z component of the angular momentum. 
This is the statement that L 2 /Q ^constant, which is 
equivalent to constant inclination in the Keplerian limit. 
In the case of B2, the result is different, but only if 73,1 7^ 
0. This leads to problems, and so in the following we 
will make this additional restriction and use the above to 
compute SQ in all cases. 

With this at hand, we can now compute E and L 
through Eq. (|260[) . To do so, every time we take a time 



derivative of the quadrupolc moment, we use the equa- 
tions of motion to reduce their order. The averaging is 
performed by taking J dt — > fdip/tp. This, of course, 
works well for the O(e ) terms (the GR contributions), 
which only depend on ip. The non-GR deformations, 
however, depend both on ip and x^ and thus, more care 
is needed. Since these are 0(e) corrections, it suffices to 
use the O(e ) piece of ip and x m the averaging of these 
terms, but at this order in e and to leading order in 1/p, 
ip = x, which then implies ip — x + ^X- After integration, 
one then obtains corrections that depend on cos 2 8x1 but 
since these quantities are slowly varying (at next order 
in 1/p, the frequencies are no longer commensurate), one 
can average over them by writing cos 2 Sx ~ 1/2. 

Let us decompose the fluxes into < E >=< E > GR 
+edE, < L >=< L > GR +eSL, and < Q >=< Q > GR 
+e5Q, where the GR fluxes are fi~5l.[ill] 

m(e) =-?4(l-e 2 ) 3/2 fl + ge 2 + ge 4 V 



\ / gr 5 p 5 y ' ' V ' 24 " ' 96 " J 

(264) 

(266) 

The correction to these fluxes are, in the B2 limit, with 
73,i = 0: 

M5E B = ~\ (1 - e 2 f 2 (712 + 2 742 ) 
b p° 



x ( l + ^W^e 4 ) , 



(267) 



7?sin6»tp , 2 x3/2 , , 

— i^T- I 1 _ e ) (^ 12 + 2 ^ 42 ) 



SL B2 = -16- 

63 2 e 4 

x H e 2 H , 

40 207 

M8Q B2 = -\^ cos 2 9 tp (l-e 2 f 2 



(268) 



x (24 + 35e 2 + e 4 ) ( 7l , 2 + 2 74 , 2 ) , (269) 



in the 6 3 limit: 



MSEb = (1 - e 2 f 2 ( 7l3 + 2743) 



, 461 2 1345 4 251 fi 

1 H e 2 H e 4 H e b 

72 288 864 



(270) 



M<5L B3 = -24^ sin0 tp (l - e 2 f /2 ( 7l3 + 2 743 ) 



361 2 7 < 
x 1 + -^e 2 + — e 



M<5Q B ; 



120 10 

— cos 2 tp (1 - e 2 ) 7l3 + 2 743 
5 p° 



(271) 
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65 2 
X,1+ 24 e 



91 



(272) 



in the B4 limit: 



208 77 

5 p s 



,3/2 



M5E Bi = (1 - ey " ( 7l4 + 2 744 ) 



1 + 


2701 


312 


96 
1248 


(1- 



2 3155 4 3779 fi 
e + — -e 4 + — — e 6 



312 



2496 

73 2 37 

— e 2 H 1 

24 96 



(273) 



M5Lb 4 



176 77 



sin 6» tp (l - e 2 ) 3/2 (714 + 2744) 



1 - 
1 

IT 



5 p 13 / 
409 



757 
352' 



!l-e 2 ) 3/2 f l + -e 2 



(274) 



M5Q B4 = -1^4 cos 2 tp (l - e 2 ) 3/2 (714 + 2 744 ) 
5 p° 



439 , 3 „ 813 4 
1 + — e 2 + — e 6 + — e 4 



104 



104 



416 



(275) 



in 



the 65 limit: 



^ 5 = -^4(l-e 2 ) 3/2 (7 15 + 
5 p y 

' 791 2 31805 4 

1 H e 2 H e 4 

.72 1728 

19277 fi 2617 8 

e 6 H e 8 

13824 



2745) 



3456 



e 2 ) 3/2 fl + ge 2 + ?Ie 4 



M5i B5 =- — -^smot 
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1 



1585 o 2493 4 

e H e 
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(276) 

(l-e 2 ) 3/2 (715 + 2745) 
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248 

^l-e^fl + Ze 2 



96 

1/2 

237 V 
496 



M<5Q Bs 



576 77 2 , 

~" fc'tp (1 - e 



— — ? cos 6» tp 
5 p' 



(277) 

2 ^ 3/2 (715 + 2745) 



, 425 2 883 4 251 B 

1 H e 2 H e 4 H e fa 

72 192 576 



(278) 



and in the CS limit: 



M^ cs = -72^ (1 - e 2 ; 



3/2 



sin tp 



1 + — e 2 + 



97 2 15065 4 1865 fi 31555 8 

" 1 e H e H e 

576 144 27648 
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73 9 37 4 N 

— e H e 
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Y 121 2 
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4585 4 
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715 2 223 4 

e 2 H e 4 
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1039 fi 
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e 2 ) 3/2 I 
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-2432 (l 



MSQ CS = -x-tth ■ n I 1 - e ) 
sm^tp 

„™„ 2. A 2347 2 1835 4 4933 fi e 8 
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4597 H 9 
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3 /2 ; „2, 



768 (l-e 2 ) 3/2 sin 2 tp 1 



(279) 



(280) 



(281) 



We can now compare these results to those obtained by 
Ryan Q [see eg. Eq. (55) in that paper]. We see that 
the CS correction to the energy flux corresponds to a 
I = 4 (hexadecapole) correction to the current multi- 
poles. Similarly, the B3 modification to the energy flux 
corresponds to a £ = 2 (quadrupole) correction to the 
mass multipoles. 

In the following section we will combine the results 



described in this section, and the precession rates derived 
in Section IIVI to obtain the leading order corrections to 
the AK waveforms in these modified spacetimes. 
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VI. MODIFIED GRAVITY AK WAVEFORMS 



A. Rate of change of eccentricity and frequency 



The AK waveform model of Barack and Cutler [12[ 
is built around the gravitational waveform from a weak- 
field, Keplcrian orbit, as computed in jTH, [T^] . The ge- 
ometry of the orbit is characterized by the eccentricity, 
e, and the radial orbital frequency, z/, while the instan- 
taneous phase of the object in the orbit is specified by 
the mean anomaly, $, which for a geodesic orbit satisfies 
d^/dt = 2ttv. The orientation of the orbit within the 
frame of reference of the BH is defined by three angles — 
the inclination of the orbital plane to the spin axis of the 
BH, A, an angle, a, describing the rotation of the orbital 
plane about the spin axis and the angle, 7, between the 
direction of pericenter and the direction L x S, where L 
is a unit vector perpendicular to the orbital plane, and 
S is a unit vector parallel to the BH spin axis. 

For a Keplcrian orbit, all of these quantities are fixed, 
but the AK model includes relativistic effects by impos- 
ing precession of the angles 7 and a and evolution of the 
orbital geometry through changes of e and v, while the 
inclination angle, A, is kept constant. This evolution is 
computed through the integration of a set of five coupled 
ordinary differential equations, which at leading order 
take the form 



d$ 
~dd 



2ins. 



(282) 



^ = ^^_ (27rM ,) 11 /3 (1 _ e2) -7/ 2 

dt 10ttA/ 2V ; V ' 



1 73 2 37 4 

1 H e 2 + — e 4 

24 96 



^ = 67rK27rMi,) 2 / 3 (l - e 2 )" 1 , 
dt 

de _ 32 7] (1-e 2 ) 3 / 2 ^ {19 t 121 
di ~~ ~ 
do 



5 M p 4 " V 6 ' 96 
?7r 2 A/^ 2 acosA(l -e 2 )~ 3/2 



To construct AK waveforms in the modified gravity 
spacetimes we have been considering here, we need only 
add corrections to these ordinary differential equations. 
The corrections to 7 and a were computed in Section HVl 
although we must change variables to rewrite these in 
terms of e and v. The corrections to e and i> are deter- 
mined by the corrections to the fluxes of energy, angular 
momentum and Carter constant computed in Section fVl 
but we must change variables to put these into the AK 
form. We will describe this below. The inclination angle, 
A, is equivalent to the inclination angle Qt p we have used 
in this paper and so the AK assumption that A = is 
consistent with keeping this angle constant, as we have 
done elsewhere in this paper. 



To convert the fluxes of energy, angular momentum 
and Carter constant into fluxes of e, p and tp we write, 
as usual, e = e K + e Se etc., where the Kerr subscript 
denotes the corresponding flux in the Kerr metric. The 
Kerr fluxes are 



Mp K 



64 (1 



(l_e 2 ) 

p3 



2\3/2 



1 



7 



32 (1-e 2 ) 3 / 2 
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#t PK = 0. 

Transforming from E, L and Q to p, e and 



dSE dSE dSE 

dp de d0 tp 

dSL dSL 8SL 

dp de dO, 

dSQ dSQ 95 
dp de 



tp ! 



(287) 

(288) 

(289) 
we find 





/ dE K dE K dE K 
dn de d9 t p 



dh v 



\ 



OLk 
de 

de 




(290) 



To leading order in p, the second matrix on the right 
hand side of this equation is 
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Rearranging Eq. (|290|) we obtain 
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and, expanding to leading order in p, M^ 1 is given by 
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As described above, for consistency with the AK model, 
we have adopted the assumption that the inclination 
angle is constant, therefore 58 tp = tpK = 0. We 
see from the above matrix that dQ K /dp/dL K /dp = 
8Q K /de/dL K /de = 2z-^/p/ sin6*t p . This result and the 
constancy of the inclination angle allows us to relate 6Q 
to SL, which gives expression (|262|) . which was used to 
derive the values of SQ presented in the previous section. 
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The constant inclination angle assumption also allows 
us to simplify Eq. (|2"9"2"|) to 
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In general, only the second term contributes to bp, while 
both terms contribute to be. 

Once bp and be have been determined, the correction 
to the rate of change of frequency can be found from 
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where we have written z> = z> K + e<5z> as usual. 

Using these expressions, we find under assumptions B2 
(and setting 73,1 = 0) 
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where the various eccentricity-dependent factors are 
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We notice that there is a problem with the e result, as 
there are terms proportional to 1 / e that do not cancel. 
This implies that circular orbits would not remain cir- 
cular under radiation reaction. It is known that circular 
orbits do remain circular in GR (42j. This proof car- 
ries over to the modified gravity spacetimes we consider 
here, and this is proven in Appendix [B] The unphys- 
ical result we have derived above has arisen from the 
assumption that the inclination of the orbital plane, 9t p , 
remains constant at leading order. It is clear that this as- 
sumption is incompatible with the physical requirement 
that circular orbits remain circular. We used the con- 
stant inclination assumption to derive Q cs from L cs for 
simplicity, and to ensure consistency with the AK model 
which also makes this assumption. For general relativity, 
a PN formula similar to Eqs. (|260[) exists to compute the 
rate of change of the Carter constant [43[ . Additionally, 
for extreme-mass-ratio inspirals, a formula for the evolu- 
tion of the Carter constant based on an expansion of the 
Teukolsky equation is also known [3, IHf . The constant 
inclination assumption is known to lead to some patholo- 
gies in numerical kludge models of extreme-mass-ratio in- 
spirals [l3| , but these arise only when this assumption is 
used in conjunction with exact Kerr geodesic trajectories, 
and the constant inclination assumption does give the 
correct leading order expression for Q. The fact that we 
have found an inconsistency only in the CS case suggests 
that the assumption of constant inclination also gives the 
correct Q at leading order for each of the other modified 
gravity spacetimes we have considered here. A full com- 
putation of Q, using either a PN expression or an expan- 
sion of the Teukolsky equation, is beyond the scope of 
the current paper, but should be explored in the future. 

In the CS case, we can derive the leading order cor- 
rection to Q in the circular orbit limit, e = 0, from the 
requirement that circular orbits remain circular under ra- 
diation reaction. That is, we can solve for the SQ CS such 
that Se C s does not contain a 1/e piece. Doing so, we find 

MSQ cs (e = 0) = -114-^7aCsin0 tp cos 2 tp , (321) 

which in fact guarantees that i5e cs (e = 0) = 0. With this 
at hand, the remaining elements become 

M6 Pcs {e = 0) = -272-i 7I?7 aCsm0 tp , (322) 
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2irM8v cs {e = 0) = -168— TyaCsin^p . (324) 

We cannot derive Q for eccentric orbits in this way, nor 
the corresponding results for Sp, Se or 59t p , although 
the eccentricity dependences of E and L given in expres- 
sion (|281j) are correct. 



B. Corrections to AK waveforms 

We now combine the above results to obtain final ex- 
pressions for the corrections to the AK waveforms. The 
final step is to rewrite the results in terms of the orbital 
frequency, is, rather than the semi-latus rectum, p. This 
can be achieved by replacing p in the preceding expres- 
sions by p K + eSp, in which 
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In this way, we get a modification to the 0(e) part of 
#7 etc., from the GR piece of these terms, evaluated at 
p K + eSp. This additional correction is sub-dominant for 
da and Sj, but can contribute at leading order to bv and 
Se. This gives final expressions for the fluxes as follows. 

Under assumptions £? 2 (and setting 73,1 = 0), these 
latter corrections arc subdominant. The final expressions 
for the corrections to the AK waveforms are 
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where g e ^ and g v ^ are as given in Eqs. (|300[) - (|302[) . 
Under assumptions Bn, we have 
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where the eccentricity-dependent factors have now been 
modified to 
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Finally, for circular orbits in the CS limit we have 

MSe cs (e = 0) = 0, (347) 
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in [28| is of Petrov type I and does not have a Carter con- 
stant. Even then, however, our modified AK waveforms 
serve as a superset of the Barack and Cutler study [28[ 
and also the Glampedakis and Babak study [37j ■ 

We note that the expressions given above depend on 
the eccentricity and inclination parameters, e and (9 tp , 
which are gauge dependent quantities. If we had con- 
structed the waveforms using a different definition of ec- 
centricity or inclination, or by parameterising the orbits 
using the energy, angular momentum and Carter con- 
stant, we would have obtained somewhat different correc- 
tions. Some care must therefore be taken when interpret- 
ing the meaning of the gauge dependent waveform pa- 
rameters p, e and 9 tp . Nonetheless, we have constructed 
the waveform corrections to ensure that quantities which 
arc measurable by an observer at infinity, specifically 
the three frequencies characterising the motion and their 
time derivatives, are consistent with the perturbed met- 
ric (modulo the omission of the "direct" contribution to 
the radiation field). Therefore, if the preceding wave- 
form corrections were now rewritten as functions of these 
three frequencies, the same, gauge-independent, wave- 
form model should result, independently of the parame- 
terisation used to derive it. 

The corrections depend on the 7^ parameters, and 
on N, which are both related to the dependence of the 
metric coefficients on the radial coordinate r. One might 
therefore expect these quantities to have a dependence 
on the choice of radial coordinate. If we were to make 
an infinitesimal coordinate transformation of the form 
r — > r + eSr(r), we could change the r dependence of 
the various metric corrections, h^ v . However, any such 
transformation would introduce hgg 7^ 0, which puts the 
metric into a different form than the one we have anal- 
ysed here. This apparent discrepancy can be understood 
from how the metric was constructed — certain choices 
were made in the form assumed for the Killing tensor 
of the perturbed spacetime, as described in Section [TlJ 
that result in the identification of the coordinate r with 
a Boyer-Lindquist-like radius. It can be interpreted as 
the circumferential radius in the equatorial plane, and 
becomes the oblate-sphcroidal radial coordinate at infin- 
ity. It is the fact that this coordinate has been fixed in 
this way that allows quantities that relate to it to be 
asymptotically measurable and to appear in the wave- 
forms. 



We can compare our results to previous results in the lit- 
erature. Barack and Cutler |28j considered a quadrupo- 
lar deformation of the Kerr metric within GR and con- 
structed the associated AK waveforms. These waveforms 
are similar to those found here in the B3 limit. In both 
cases, 5u and £7 scale as v 5 and v 1 ^ [see Eqs. (6) and 
(7) in [28|]. Differences arise in the evolution equation for 
Se and 5a, which are probably due to the fact that the 
metrics used here and in [28j cannot be mapped into each 
other. This is because our metric is constructed such that 
an approximate Carter constant remains, while that used 



C. Example Waveform 

To illustrate the modified gravity waveforms we show 
an example in Figure [TJ This plot shows one of the two 
low-frequency Michelson res pon se channels of the LISA 
detector, hi, as described in pjj], for an EMRI occurring 
in the Kerr metric, and for EMRIs occurring in a de- 
formed metric with a $2-type deviation. We use two dif- 
ferent values for the size of the deformation of the metric, 
£(71,2 + 274,2) = 0.01 and 6(71,2 + 274,2) = 0.1. The other 
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waveform parameters were taken to be M = 10 6 Af Q , = 
10M Q , S/M 2 = 0.7, e = 0.25, A = tt/4, u = 0.00177, 
6 S = tt/4, </>5 = tt/2-34, fljf = tt/8, D = lGpc and 
4>k = otQ = 70 = $o = Oj using the notation of the an- 
alytic kludge model [jjjj. We show two sections of the 
waveform, each of duration 10 4 s, taken at the beginning 
and the end of a 2 week long observation. 

We see that the B2 waveform for the smaller perturba- 
tion is almost indistinguishable from the Kerr waveform 
for the first 10 4 s, while for the larger perturbation the 
waveform can be seen to be beginning to drift out of 
phase already by the end of the first segment. All three 
waveforms arc clearly out of phase in the second segment 
and would therefore be clearly distinguishable in a LISA 
observation of two weeks duration. We note, however, 
that we have fixed the parameters for all three waveforms 
to the same values, and have not made any attempt to 
improve the overlap between the modified gravity wave- 
forms and the Kerr waveform templates by making small 
adjustments to the parameters. Such parameter adjust- 
ment would better approximate the process of parameter 
estimation through matched filtering that will be used to 
analyse LISA data in practice and it is likely that this 
would significantly improve the phase coherence between 
the Kerr and modified gravity waveforms. Nonetheless, 
this example illustrates that LISA should be able to place 
some kind of constraint on the size of any deviations from 
GR of this form that would be consistent with EMRIs ob- 
served by LISA. We leave a full study of such constraints, 
accounting for parameter correlations, to the future. 



VII. DISCUSSION 

We have here taken the first steps toward the con- 
struction of ppE waveforms for EMRIs. We have taken 
the parametrically deformed, bumpy metric for modi- 
fied gravity theories described in [34j and calculated the 
geodesic equations in second and first-order form, the 
orbital frequencies and the implicit deformations to the 
fluxes of energy, angular momentum and Carter constant. 
With all these ingredients, we then explicitly described 
how to construct ppE AK waveforms. 

A natural follow-up to this work would be to study 
the accuracy with which space-borne detectors, such as 
LISA, could constrain such ppE deformations. It was 
in the context of such a study for EMRIs in a Kerr 
background that the original AK model was developed 
by Barack and Cutler [12|]. Now that we have obtained 
the leading order corrections to the AK waveforms, we 
can follow a similar analysis to ask with what precision 
LISA or a LISA-like mission might be able to measure 
the bump parameters Jn,m that characterize the devia- 
tions from GR. The modifications described here could 
also be incorporated into the numerical kludge waveform 
model [lH for a similar study. Higher-order PN correc- 
tions have already been included in that framework for 
Kerr EMRIs. 
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FIG. 1: Example of a modified gravity waveform. In each 
panel, we show the hi Michelson response for EMRIs in the 
Kerr metric (red solid line) and in a £?2-type modified gravity 
deformed metric with £(71,2+274,2) = 0.01 (green dashed line) 
and £(71,2 + 274,2) = 0.1 (blue dotted line). The two panels 
show the first (upper) and last (lower) 10 4 s of a 2 week long 
waveform. 



Another possible application of these results would 
be to study how these GR deformations could be con- 
strained with observations in the electromagnetic spec- 
trum, such as from the emission spectra of AGN or iron 
line profiles |46| . The emission properties of discs are 
determined in part by the geodesic structure of the BH 
spacetime, which is characterized by the orbital frequen- 
cies and precession rates computed in this paper. Wc 
note, however, that such an analysis would have to ac- 
count for the gauge dependence of these results, and work 
with gauge- invariant quantities, for example by rewriting 
the expressions in terms of the frequencies of the motion, 
or some other asymptotic observables. We have not wor- 
ried about this here, as the goal was to provide a frame- 
work in which to compute waveform observables and per- 
form a data analysis study. The specific paramctrisa- 
tion used to describe the waveform family does not affect 
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those aims, although care will have to be taken when 
interpreting the observations that are made using these 
template waveforms. 

The construction of ppE waveforms for modified grav- 
ity bumpy metrics is not yet complete, as we have here 
neglected explicit or direct corrections to the radiation- 
reaction force. As shown in [3l|, corrections to the con- 
servative and dissipative sectors of any theory lead to 
waveform modifications that might be degenerate. Al- 
though this is definitely the case for comparable mass, 
quasi-circular inspirals, the generality of EMRI orbits 
might break these degeneracies. Thus, the work done here 
can be considered as a solid first-step toward generic ppE 
EMRI waveforms, where future studies should concen- 
trate on the addition of dissipative effects. 

Complete results for the corrections to the radiation- 
reaction terms in CS modified gravity were only pre- 
sented for circular orbits, as we had made the simplifying 
assumption that the change of the orbital inclination due 
to radiation reaction was a sub-dominant effect. This 
assumption made the computation of the corrections to 
Q more tractable, but it appears not to be valid in CS 
modified spacctimes. We were able to present results for 
circular orbits by imposing the physical requirement that 
circular orbits remain circular under radiation reaction. 
An extension of the Q calculation to the case of generic 
orbits should be pursued in the future. 
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Appendix A: Useful integrals 

Computation of the frequencies has relied in various 
places on the evaluation of integrals of the general form 



/„ = / (1 + e cos ■0)"d'0 
Jo 

and we describe here how these may be computed. 
For n > we can expand the integrand as 
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in which [a^J denotes the largest integer smaller than x. 
On integration over the interval [0,7r], all the terms in 
the summation give zero. The second term is only non- 
zero when k is even and therefore only the even powers 
of cos?/' contribute to the result. We therefore find 
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in which n\ = n(n — l)...l is the usual factorial. If n < 0, 
the integral can be evaluated by first writing cos ip = 
cos 2 (-0/2) — sin 2 ( , i/'/2) and then using a substitution t — 
tan(?/)/2). This reduces the integral to 
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Finally, writing (l + e)/(l — e) = £we can see through 
integration by parts 
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The integrals we will need for the various special cases 
we discuss in this paper are 

7T 7T 



~* (!_ e 2)3/2' 
I = h = 7T, h = 7T ( 1 



(l_ e 2)l/2' 
2 * 



h = 7T 1 



3_! 
2 



h = 7T 1 + 3e 



3e 4 



(l + ecosi/;)™ = ^('^ N )e fe cos' £ V- (A2) 



I 5 = 7T 1 + 5e 



15e 4 



(A6) 
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Appendix B: Circular orbits Remain Circular 

In thi s ap pendix, we study the proof of Kennefick and 
Ori [42|, |47[ that circular orbits remain circular, but ex- 
tended to the non-Kerr backgrounds analyzed in this pa- 
per. Let us begin by recasting the radial equation of 
motion as an equation for the carter constant 



Q = H GR (r, E, L) + SH{r, E, L) - Ai 



(Bl) 



where as usual A = r 2 — 2AIr+M 2 a 2 , while u r = g rr u r 
p 2 r/A. The potentials H GR and SH are defined via 

H GR = A' 1 [E(r 2 + M 2 a 2 ) - M 2 aL] 2 
- M 2 (L - aE) 2 - r 2 , 
SR 



SH = 



AM 2 ' 

where the SH in the CS limit becomes 
5 MEL (MX 



(B2) 



(B3) 



The action of an external force on a particle will force 
the evolution of the Carter constant. When this force is 
the radiation reaction force, then 



Du a 

~d7' 



(B4) 



where D stands for covariant proper-time differentiation. 
The Carter constant then evolves as 



4-g-F. 
ou a 



(B5) 



Using the instantaneous circularity condition, u r = 0, 
plus the chain rule, the above equation becomes 



Q = Q,eE + Q, l L. 



(B6) 



Using Eq. (|B1[) , the evolution of Q for a circular orbit is 
simply 

Q = H G £E + H% R L + (SR E E + 5R >L L) , (B7) 

which in the CS limit becomes 



rB • 5 M 2 EL ( M\ 2 ( E L 
Q = H%»E + HTL +l —ad-) l- + J - 



(B8) 



Let us now compute the rate of change of the Carter 
constant required to take a circular orbit into a new cir- 
cular orbit. Consider then Eq. (|Blj) . but solve for Au 2 . 
This new equation will then resemble that of a parti- 
cle in an effective potential, given by W(r, E, L, Q) = 
H Gn + 5H — Q. The stable circularity condition is equiv- 
alent to requiring that W = and W. r = 0. 

If the radiation reaction force acts on this instanta- 
neously circular geodesic for a short time St only. After 
this, the orbit will change to a new geodesic with pa- 
rameters X' = (r' , E' , L' ,Q'). Denoting the change by 
ST = X — X', the change in the potential is then 



SW = W.eSE + W. L 5L + W.qSQ , 



(B9) 



where we have used the fact that W <r 
orbits. 



for circular 



For this equation to hold before and after the action 
of the radiation reaction force, we must have SW = 0, 
which allows us to solve for SQ. After dividing by St and 
taking the St — > limit, we find 



Qc 



We ■ 
W.n 



Wj 
Wc 



-L. 



(BIO) 



Using the definition W, and assuming that W t Q 
wc then find 



Q 



_ rrGR 771 
circ — E & 



1 



SR.eE + sr l )l 



AM 2 

(Bll) 

The assumption that W t Q = — 1 relies on the condition 
that SR be independent of Q. This condition is satisfied 
in CS gravity and in all Bn limits. It is not, however, 
satisfied if one modifies the £>3 limit by letting 71.2 7^ 0, or 
if one modifies the B4 limit by letting 7^3 ^ 271.2 and if 
one modifies the B$ limit by letting a 2 7i,2 7^ 27^3 — 71,4. 



The fact that Q C irc in Eq. (|B11[) is identically equal 
to Q in Eq. (|B7I) means that even after the radiation- 
reaction force acts on an instantaneously circular orbit, 
the induced, instantaneous evolution of the orbit's Carter 
constant is such as to maintain instantaneous circularity. 
This result is identical to what is found for geodesies in a 
Kerr background. We can trace back the reason for this 
to the fact that the potential SH (or equivalently SR) 
is independent of the Carter constant in all the cases 
studied here, including the CS one. 
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